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ABSTRACT 


The work presented here is a study of the behavior of the 
electronic impedance under transient conditions for a type of Impatt diode. 
This study includes both analytical and experimental investigations and 


is applicable to Read model and one-sided abrupt junction Impatt diodes. 


An improved analytical model of the diode has been formulated 
and evaluated. This model allows the inclusion of the space charge 
effect associated with the avalanche region as well as individual ionization 
rates and drift velocities for both types of carriers. The model spans 
the range of applicability of both small-and large-signal analysis and 


can be used for steady-state or transient operating conditions. 


The time dependence of the Impatt diode's junction temperature 
is measured to within five nanoseconds of the onset of the bias pulse and 
a model for the transient heat flow in the diode's semiconductor junction 
is developed. This heat flow model gives close agreement with the 
experimental results and is applicable for a range of pulsed operating 
conditions. Also, a means of incorporating changes in the junction 


temperature into the analytical model is evaluated. 


A method for the measurement of the diode's electronic impedance 
at microwave frequencies and under transient operating conditions is 
developed. This method employs the principle of reinserting the carrier 
frequency to measure the complex reflection coefficient as a function 
of time and is referred to here, as the Carrier Reinsertion Method. With 


this technique, the diode's electronic impedance is measured during the 
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buildup and decay of bias pulses having risetimes of five nanoseconds. 
There is good agreement between these measured results and those 


provided by the analytical model. 


The analytical model of the Impatt diode and the results of the 
transient electronic impedance measurements are utilized to construct a 
computer model of the diode. This model indicates that the transient 
characteristics of the diode's electronic impedance have a relatively 
short duration and are largely confined to the real part of the electronic 
impedance. Also indicated is the dependence of the transient effects 
upon the rise times of the bias voltage pulse and on the frequency and 
magnitude of the R.F. voltage impressed across the diode. The significance 
of these transient effects with respect to Impatt diodes used in R.F. 


amplifiers and oscillators is discussed. 
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CHAPTER I 
INTRODUCTION 


In recent years, the use of semiconductor devices in micro- 
wave communications has increased significantly. Within the past 
several years, solid-state microwave oscillators and amplifiers have 
demonstrated economic and/or performance advantages over microwave tube 


esd : ; 2 ; 
(1, ) and as a result are being incorporated into microwave 


components 
systems. Two important semiconductor devices used for the generation 
and amplification of microwave frequencies are the Transferred Electron 
and Avalanche diodes. The Transferred Electron devices are most suited 
for low-noise, low-power applications, while the Avalanche devices, at 
the present time, are the most powerful solid-state source of micro- 


(3) 


wave energy for continuous wave (CW) applications‘~’. 


It is unlikely that the avalanche diode,or any other solid- 
state device, will simply replace tube devices in the near future. 
Rather, solid-state technology will be applied selectively to perform 
specific system functions which freely capitalize on the advantages of an 
all solid-state approach. For this reason and to ensure successful 
implementation in the new systems concepts, a precise knowledge of 
the characteristics and limitations of these solid-state devices is 


desirable. 


This research project is concerned with studying the transient 
behavior of pulsed Avalanche diodes in their principal mode of operation. 


Avalanche diodes may operate in several modes \°?, some of which are 
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referred to as TRAPATT, parametric and Thermal, with the principal 

mode being referred to as the Impatt mode (Impact ionisation Avalanche 
and Transit Time). Mathematical models describing the steady- 

state behaviour of this diode in different modes of operation have been 
publ ished (42526) , but very little work has been carried out, so far, 


on the transient behavior, particularly during pulsed operation. 


The buildup of oscillations in an Avalanche diode oscillator 
and the operation of Avalanche diode amplifiers under varying signal 
conditions, are actually a study of the behavior of the diode's 
electronic impedance under transient conditions. The transient 
conditions may be those of initial turn on-off or the rapid-step changes 
in the level of R.F. power. The term "electronic impedance" refers 
Only to the impedance of the active part of the semiconductor, while 
the term "diode impedance" refers to the combined effect of the 
electronic impedance, the impedance of the passive part of the semi- 
conductor and the diode package impedance where applicable. This 
electronic impedance of an Impatt diode is dependent upon the d.c. 
bias voltage, as well as the frequency and magnitude of the R.F. voltage 
across the diode. Consequently, some special problems arise in the 


design of Impatt diode amplifiers and oscillators. 


Solid-state circdsts! @ capable of small-signal amplification 
at microwave frequencies, such as tunnel diode - and parametric 
amplifiers, have been in existence for many years. It has been within 


the last decade that semiconductor devices with potential high power 
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amplifying properties, such as Avalanche diodes, have become available. 
Several interesting mireut concepts have already been realized using 
Impatt diodes that promise to open up new areas of system 

applications (7), Due to the negative resistance mechanism of 

the Impatt diode both the real and imaginary parts of the diode's 
impedance are strongly dependent on the d.c. bias current and R.F. 
signal levels. Shifts in the frequency and magnitude of the maximum 


(8) 


gain occur for different incident R.F. levels, which leads to some 
distortion in the amplification and injection locking characteristics 

of reflection amplifiers. A more thorough understanding of the transient 
behaviour of the diode's electronic impedance would be helpful in 


designing these amplifiers. 


(9) 


In the communications industry there is a growing need 

for an inexpensive, short range, common carrier facility which has a 

high message-carrying capacity. This system would serve such purposes 

as intercity computer - to - computer links where the cost of laying 
cable would be prohibitively high. A two-level pulse code modulation 
(pcm) microwave repeater system with a bias modulated Impatt diode as 

a transmitter stage would fulfill this need. A bias modulated oscillator 
combines the oscillator and modulation functions into one unit. This 
system has significant cost advantages as well as requiring a minimum 

of hardware. Two ways in which two-level PCM can be produced are 


frequency modulation (FM), with a self-deviating oscillator, and amplitude 


modulation (AM) of the R.F. carrier. 


In the self-deviating Impatt oscillator, the R.F. carrier is 


continuously produced at a constant amplitude and the bias modulation 
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Changes the electronic impedance of the diode and hence the frequency 

of the R.F. oscillations. When the complex electronic impedance of the 
diode changes, both the real and imaginary parts change, which con- 
tributes to both amplitude and frequency changes in the R.F. oscillations. 
Consequently, some distortion from amplitude-frequency conversion can 

be expected. This, coupled with other characteristics of the FM 
modulation of Impatt diodes that are not fully Micerereod: ie 
necessitates a better understanding of the Impatt diode's electronic 
impedance under transient conditions, particularly at the higher 
modulating frequencies. 


(10) ys 


For amplitude modulation of the R.F. carrier 
procedure is commonly referred to as Amplitude Shift Keying (ASK), where 
a binary system of either on or off can be used for the oscillator. For 
this ASK system, the R.F. fields have to grow ina cavity,and this 
buildup will depend on the circuit parameters as well as on the electronic 
impedance of the diode. The diode's electronic impedance is a function 
of the modulating bias voltage and also the magnitude of the R.F. in 
the cavity. Therefore, a knowledge of the behaviour of the electronic 
impedance under transient conditions is necessary for proper circuit 
design. Further, the pulse lengths may be appreciable and the effects 
of the thermal properties of the diode junction on the performance of 
the oscillator will become important. Since the diode is pulse operated, 
the instantaneous junction temperature varies with time. This variation 
in temperature will cause both amplitude and frequency modulation 
CanaC terres a that may be detrimental to the system operation as 


a whole. Therefore, the transient thermal characteristics of the 
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diode's junction temperature should be known for design purposes. 


The discussion of Impatt diodes is usually separated into 
small-and large-signal analysis. Small-signal results predict the range 
of current and frequency in which an Impatt diode will oscillate if 
placed in a proper resonate circuit. These results also indicate the 
magnitude of the maximum negative impedance which can be generated 
by the diode. This information is very useful during the initial 
stages of device evaluation for oscillator and amplifier circuits. 
However, small-signal analysis does not provide sufficient information  !2? 
concerning maximum power output, efficiency, saturation and stability 
characteristics which are important in large-signal oscillator and 
amplifier designs. Thus, the requirement for and the name, large- 
Signal analysis. The analytic solutions for steady-state small-signal 
analysis are available in closed form due to the many approximations 
that can be made. While simplified analytic solutions for the steady- 
state large-signal case may be presented in closed form, the usual 
procedure is to obtain computer solutions to less simplified 
describing equations. A study of the transient behaviour of the elect- 
ronic impedance of an Impatt diode will span the range of applicability 
of both small - and large-signal analysis, and incorporate the time- 
dependence of the d.c. bias voltage pulses, with the most probable 


form of solution being a computer model of the diode. 


The presentation of this work is divided into several parts. 
Chapter II discusses the general theory of operation and the particular 


device structure chosen to represent the Impatt diode. Chapter III 
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describes the mechanisms and the assumptions made in the physical 

model along with the fundamental equations and boundary conditions 

that mathematically determine both the steady-state and transient 
problems. Chapter IV describes the details of the analytical formulation 
and computational techniques used in the solution of the time dependent 
equations for the Impatt diode. Chapter V presents the analytical and 
experimental investigation of the transient heat flow in the diode's 
junction. Chapter VI covers the experimental techniques developed for 
the measurement of the electronic impedance of the Impatt diode 

under transient conditions. Chapter VII provides a comparison of the 
computational model and the experimental results. Chapter VIII 
presents a discussion and conclusions drawn from the results of this 


study. 
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CHAPTER II 
IMPATT DIODES 


When a semiconductor p-n junction is reverse-biased, a 
Space charge region is formed that is depleted of mobile carriers and 
only a small reverse saturation current flows. When the reverse bias 
voltage is increased until it exceeds a critical value, the junction 
breaks down and a large current will then flow for only a small increase 
in bias voltage. This large current is the result of avalanche 
breakdown, when internal secondary emmission causes a multiplication 
of electrons and holes in the depletion region of the junction, 
The critical level of reverse bias voltage, at which the multiplication 
begins ,is referred to as the breakdown voltage. The p-n junction, in 
the avalanche breakdown condition, may exhibit a negative-resistance 
characteristic in the microwave frequency range. A diode (p-n 
junction) with negative-resistance terminal characteristics, has the 
capability of converting d.c. energy to R.F. energy and thus an R.F. 
voltage present across the diode will grow in amplitude until limited 


by circuit and diode losses. 


Avalanche diodes can be designed and operated in various 
modes. The first and most important mode ''3) is the IMPATT mode. 
(Impact ionization Avalanche and Transit Time). This mode produces its 
negative-resistance terminal characteristics by having the carrier 
density lag the applied electric field by more than 90°, which is also 
the phase difference between the fundamental R.F. component of the 


diode's external current and the R.F. voltage across the diode. 
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2.1 Basic Principles of Operation 


The more common materials used in the fabrication of avalanche 
diodes are Silicon, Germanium and Gallium Arsenide. Due to the basic 
physical phenomena involved in the generation of the negative-resistance 
in this device, almost any semiconductor material can, in principle, be 
employed in the construction of avalanche diodes. The Impatt mode 
of operation involves two basic properties of carrier transport in 
solids at high electric fields, namely those of avalanche multiplication 


and transit-time delay. 


The breakdown of the p-n junction is caused by an avalanche 
multiplication process involving both electrons and holes (as opposed 
to other processes such as zener breakdown and tunneling). When the 
electric field becomes sufficiently high, a charge carrier (electron 
or hole) can acquire sufficient energy from the electric field, between 
collisions, to dislodge a bound electron into the conduction band, 
during a collision, thus also creating a hole in the valence band. 

This action is referred to as impact ionization and occurs for magnitudes 
of the electric field greater than about 10° V/cm (for Silicon). The 
probability of this electron-hole pair creation, or ionization rate, 

is a sensitive function of the electric field strength. The ionization 
rate can increase as much as five orders of magnitude for an increase 

in the electric field strength of only two or three times (4) When 

the electric field changes periodically with time around an average 
value, the generation rate of carriers follows the field change almost 
instantaneously. However, because the total generation of carriers 


depends upon the number of carriers present, the total number of carriers 
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generated does not change in unison with the electric field and keeps 
On increasing after the electric field has passed its peak value. The 
total number of carriers peaks and starts to decrease when the field 
has decreased from the peak to the average value. Thus, bunches of 
carriers are formed and the a.c. variation in the number of carriers 


lags the a.c. component of the electric field by 90°. 


When the electric field is low, the average drift velocity 
of the carriers is directly proportional to the electric field and 
the proportionality constant is called the low field mobility. As 
the field strength increases, carriers interact more strongly with 
the lattice so that the average drift velocity falls much below that 
of a linear projection from the low field values. For electric field 
strengths above about 10° V/cm(for Silicon) sthe drift velocity 
approaches a limiting value and becomes independent of the electric 
field. This limiting value is referred to as the scattering-1limited 
drift velocity. When a bunch of carriers is injected into a region 
where the electric field is such that avalanche multiplication does 
not occur and yet the field is sufficiently high that scattering 
limited drift velocities apply, the transit-time of the bunch will be 
directly proportional to the width of the region. This transit-time 
can be used to delay the collection of the bunch of carriers and thus 
introduce a phase shift between the fundamental R.F. component of the 
external current, due to the transiting bunch of carriers, and the 


R.F. voltage across the transit-time region. 
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The above discussion of avalanche multiplication and 
transit-time delay has been for ideal conditions. In practice, the 
effects of space charge and non-scattering-limited drift velocities 
must be taken into account. In the region where avalanche multi- 
plication occurs, the bunch of carriers causes the electric field to 
be depressed. This depression of the electric field due to the space 
charge effect of the carriers, results in the field strength falling 
below the average value before it would normally do so, when affected 
by the R.F. voltage alone. Therefore, the phase lag of the bunch of 
carriers, due to the avalanche process, is decreased and will be less 
than the optimum 90° with respect to the R.F. voltage. This space 
charge effect will increase with increasing magnitude of the bunch of 
carriers. When the electric field in the transit-time region falls 
below that required to maintain scattering=limited values of drift 
velocity, the drift velocity becomes dependent on the magnitude of the 
electric field. When this occurs, dielectric relaxation tends to 
smooth out the carrier bunching. This "debunching" reduces the magnitude 
of the current external to the transit-time region and hence reduces 


the power output. 


Generally, Impatt diodes consist of an avalanche region, 
where carrier multiplication occurs, appended to a transit-time (drift) 
region where carrier multiplication does not occur and the drift 
velocity of the carriers is assumed to be scattering limited. In 
the avalanche region holes and electrons are produced in practically 
equal numbers and since the multiplication is very sensitive to the 


magnitude of the electric field, it is here that the carrier transport 
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current is modulated by the R.F. field. The contribution of the 
avalanche process to the phase difference between the diode's external 
conduction current and the R.F. voltage across the diode will be 
something less than 90° due to the space charge effects of the bunches 
of carriers. As these bunches of carriers "drift" across the transit- 
time region, the phase difference added to that of the avalanche 
process 1s approximately equal to one-half the transit-time. The 
working region of the diode consists of both the avalanche and drift 
regions, with the total delay in the conduction current being the sum 


of the avalanche and transit-time processes. 


2.2 Types of Impatt Diodes 


The basic members of the Impatt diode family 4 are the Read 
diode, one-sided abrupt p-n junction, linearly-graded p-n junction and 
the p-i-n diode. Only the Read diode and the one-sided abrupt p-n junction 
are of interest for this work. The Read diode is of interest because 
it is a convenient structure for analytical purposes and contributes 
to an understanding of the dynamic operating characteristics of the 
Impatt diode. The one-sided abrupt p-n junction is of interest because 
it is the type of Impatt diode that is most commonly manufactured and 
also is the type which will be used in the experimental portion of this 
work. The linearly-graded p-n junction and the p-i-n diode are not of 
direct interest for this work and have only been included here for 


completeness. 


Fig. 2.1 shows the doping profile, electric field distribution 


and ionization integrand, at the breakdown condition, for an idealized 
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Read diode (pn v n or its dual n° 8) Bey The corresponding a.c. 
voltage and current waveforms are shown in Fig. 2.2. For the Read 
model, all of the avalanche process is confined to the n region of 

the pn v n diode (in the dual model - the p region) where the actual 
avalanche region is defined by the ionization integrand. It is 

assumed that this avalanche region is sufficiently narrow that its 
transit time is a negligible portion of an R.F. cycle. The diode is 
designed such that. the electric field in the drift region (the 

intrinsic region) is sufficiently high to ensure that the carriers 
travel at their scattering limited drift velocities. In Fig. 2.2, 

the idealized a.c. operation of the diode is depicted. The a.c. 
component of the electric field (ey oot modulates the ionization 

rates for the carriers to forma pulse of carriers in the avalanche 
region. As this pulse of carriers is formed the space charge effect of 
the carriers (AE) also increases and, in practice, the peak of the pulse 
of carriers will occur before the electric field has returned from its 
peak to its average value. Thus, the phase shift in the avalanche region 
will be somewhat less than the 90° shown in Fig. 2.2. After the pulse of 
carriers has been formed, it travels across the transit-time region at 
the scattering limited drift velocity. For the duration of this transit 
time, of the pulse of carriers, a steady current flows in the external 


circuit (I_.). Once the pulse of carriers reaches the collecting 


eX 
terminal of the drift region the electric field in the avalanche region 
rises above the average value and the whole process repeats itself. 
As shown in Fig. 2.2, the fundamental R.F. component of the external 


current is delayed with respect to the R.F. voltage by the phase lag 
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FIG. 2.2 VOLTAGE AND CURRENT WAVEFORMS FOR IMPATT DIODE 
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inherent in the avalanche process plus the phase lag of one-half the 
transit-time of the drift region. In this manner the negative- 


resistance terminal characteristics of the Read diode are formed. 


In Fig. 2.3, the doping profile, electric field distribution 
and ionization integrand for the breakdown condition of an idealized 
One-sided abrupt p-n junction are given. The gradient of the electric 
field is sufficiently steep that the ionization integrand and hence 
the width of the avalanche region are localized near one end of the 
depletion region. Therefore, the total depletion region can be 
divided into separate avalanche and drift regions, with the analysis 
being carried out in a similar manner to that of the Read diode. A 
prominent distinction that occurs between these two types of Impatt 
diode is in the length of the drift region. For the field in the 
drift region to be sufficiently high to ensure that the carriers travel 
at their scattering limited velocities, the drift region may not be 
long enough to give a full 90° phase lag (at the microwave frequency in 
use) as in the Read model. Also, any portion of the drift region 
where the drift velocity falls below its scattering-limited value will 
contribute a debunching effect to the pulse of carriers, with a 
corresponding decrease in the magnitude of the external current. There- 
fore, the one-sided abrupt p-n junction will be less efficient than 
the Read diode. However, the one-sided abrupt p-n junction is more 
readily manufactured and therefore has a distinct practical advantage 


over its Read model counterpart. 
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(c) 
FIG. 2.3 ONE-SIDED ABRUPT P-N JUNCTION 
DOPING PROFILE (ptn) NET IMPURITY CONCENTRATION 
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2.3 Fabrication of Impatt Diodes (15) 


As mentioned previously in Section 2.2, the more common 
materials for the fabrication of Impatt diodes are Silicon, Germanium 
and Gallium Arsenide. The material parameters such as ionization rate, 
dielectric constant and thermal conductivity will alter the 
characteristics of the Impatt diode when constructed from different 
materials. Germanium has the highest ionization rate and therefore the 
lowest breakdown field. This is an advantage because it reduces the 
necessary input power for Impatt operation. The ionization rates for 
Gallium Arsenide and Silicon are about the same. Germanium has the 
largest dielectric constant and will have lower impedance levels and 
also higher bias currents, so this is a disadvantage. Silicon and 
Gallium Arsenide have similar dielectric constants. The thermal 
conductivity of Germanium is only about 30% that of Silicon, with 
Gallium Arsenide having a lower value of thermal conductivity than 
Germanium. Because of the relatively high values of reverse bias 
voltage required for Impatt diode operation, an appreciable amount of 
heat is generated in the p-n junction and a high value of thermal 
conductivity is advantageous. The technology for material preparation 
and processing is more highly developed for Silicon and this, along with 
its high thermal conductivity, make Silicon the most commonly used 


material for the fabrication of Impatt diodes. 


Fabricating a Read model of the Impatt diode is more difficult 
than fabricating a one-sided abrupt p-n junction model. In the 


+ 
formation of the Read diode, using a pny n as an example, a double 
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diffusion technique is required where a low concentration of donors 

is diffused deeply and then a shallow aifFosi oh of a high concentration 
of acceptors is carried out. This is not the only technique for the 
fabrication of Read diodes but is a typical illustration of the 
complexities involved. In the fabrication of abrupt p-n junctions 

only a single diffusion is necessary,namely that of a shallow diffusion 
of a high concentration of acceptors. However, care must be taken 

that the junction depth ( n region) is less than the space-charge 
(depletion) region for the breakdown voltage associated with the doping 
concentration of the n region. Otherwise, part of the high resistivity 
n region will not be swept clear of mobile carriers when the reverse 
bias voltage has reached the breakdown level and also will have a 
magnitude of electric field that will be too low to maintain scattering 
limited drift velocities. In theory, the Read diode should be more 
efficient than the one-sided abrupt p-n junction however, according 

to the results reported in the exscerauiaa, after fabrication, the 


abrupt junction diodes may outperform their Read diode counterparts. 


In avalanche diodes, the d.c. power which is not converted 
to microwave power, generates heat in the high field region near the 
junction. As the junction temperature increases, the breakdown 
voltage increases and the ionization rates decrease. Under conditions 
of uniform current density, the center of the junction area will be 
hottest and the temperature dependence of the ionization rates and 
breakdown voltage will tend to decrease the current density in the 
diode'scentral area. Thus, the temperature distribution across the 


diode junction is more uniform but the current distribution will be non 
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uniform and tend to concentrate near the outer regions of the junction 
area. It is therefore advantageous to design the diode with the 
junction near the semiconductor surface and to mount the diode with 
this semiconductor surface in contact with a low thermal resistance 


path to the heat sink. 
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CHAPTER ITI 
PHYSICAL AND MATHEMATICAL MODEL 


3.1 Physical Model 


The work reported here is concerned with the behaviour of 
Impatt diodes. The characteristics of the semiconductor material and 
the formation of the p-n junction are of interest in so far as they affect 
the operation of the Impatt diode. Some of the physical characteristics 
of semiconductors that may be of interest and the assumptions made, will 
now be discussed. 


By ole carrier Peconbi nation ee 


Under conditions of thermal equilibrium, carrier generation is 
balanced by carrier recombination. When thermal equilibrium is disturbed, 
the mechanisms that restore equilibrium are the drift and diffusion of 
carriers into and out of the region, along with recombination. The lower 


(14) 


limit of the lifetime of minority carriers in silicon is of the order of 


1078. seconds. The transit time of the carriers in Impatt diodes is of the 


=10 seconds, for diodes designed to function below 10GHz. There- 


order of 10 
fore, although the Impatt diodes may operate with high current densities, 
recombination should be of small importance and will be assumed to be 
negligible for purposes of this study. 


3.1.2 Surface Boundary condi ene we 


In this analysis, it is assumed that thermal equilibrium prevails 
at the contacts, so that the carrier density at the contacts is constant 


and equal to its thermal equilibrium value at all voltages of interest. 
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Along with the condition for thermal equilibrium, it is also assumed that 
charge neutrality prevails at the contacts. These two assumptions are 
generally regarded as the basis for having ohmic contacts. Some recombination 
may also occur at the surfaces but this is assumed to be negligible in 


this study. 


Sibi bias! Degeneracy!'®) 


For purposes of this study, Boltzman's transport equation and 
Boltzman statistics are assumed to be valid. These assumptions simplify 
the mathematical analysis as changes in transport phenomena can be derived 
from Boltzman's transport equation and the validity of Einstein's relation 


(18) the difference between Fermi 


can be assumed. As discussed by Shockley 
and Boltzman's statistics is negligible when the predicted carrier 
concentrations are small compared to the effective number of states in 
the conduction and valence bands. For Silicon, Boltzman statistics are 
valid up to carrier densities of 10! 8¢m73,, (16) which is much higher than the 
normal doping levels in Impatt diodes. Therefore, the assumption of a 
non-degenerate semiconductor will be made. 


3.1:4 Carrier Mobility 6.19 


The mobility of electrons and holes is determined by the type 
of dominant scattering. The two most significant forms of scattering are 
those caused by ionized impurities and those caused by mechanical vibrations 
of the crystal. Scattering by ionized impurities is coulomb scattering and 
is most effective at low velocities, while mechanical vibrations cause 


acoustic and optical phonon scattering. For high electric fields (~10°v/em) , 
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which is the normal operating range for Impatt diodes, the carriers 

gain appreciable energy from the electric field. As the field increases, 
the average energy of the carriers also increases, and they acquire an 
effective temperature which is higher than the lattice temperature. 

These high electric fields heat the carriers at a sufficient rate that 
the mobility decreases in proportion to the rate of increase in the 
electric field, causing the carrier drift velocity to become scattering 


limited. 


GAinDe MI Cro 1asmas (1622021) 


As a result of defects in the crystal structure of the semi- 
conductor, small areas are formed that behave as a gaseous plasma. These 
Phenomena are often referred to as microplasmas. When the Impatt diodes are 
reverse-biased, these microplasma areas reach the avalanche condition at 
an applied voltage that is a few volts lower than the remainder of the 
junction. The resulting avalanche breakdown starts at several locations 
in the junction and gradually spreads, with increasing voltage, until the 
whole junction avalanches. This gradual and somewhat unpredictable start 
of avalanche breakdown can be avoided by biasing the Impatt diode, such 


that a steady d.c. current of one milliampere is maintained in the diode. 


3.2 Validity of One-Dimensional Analysis 


In the use of more than one-dimensional analysis, all the physical 
parameters for a particular diode would be required and the complexity 
of the describing equations would be greatly increased. The Impatt diodes 
available for this work are not adequately characterized to warrant the 
greatly increased effort and expense required to carry out an analysis 


in more than one-dimension. 
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Further, the physical construction of the Impatt diodes tends 
to accentuate one dimensional behaviour. In a typical, low power, Impatt 
diode, the cross-sectional area is about 107" em? and the length of the 
activity region is 3-5 microns (for oscillations around 7 GHz). This 
results in a diameter to length ratio of around 30:1. Since it is assumed 
that avalanche multiplication occurs uniformly over the cross-sectional 
area, the carriers will travel mainly directly through the diode. Only 
under very high current densities, where carrier mutual repulsion and non- 
uniformity of temperature become appreciable, will the lateral spread of 
carriers be significant. The Impatt diode may be operated at high current 
densities and some inaccuracy will be inherent in the use of one-dimensional 


analysis, however this should be a small effect. Therefore, one-dimensional 


analysis will be used throughout this work. 
3.3 Mathematical Model 


Having made the previous simplifying assumptions for the physical 
processes occurring inside the device, the next step is to obtain math- 
ematical equations which will model the physical behaviour of the device. 


These mathematical equations will be written in one-dimensional form. 


3.3.1 Fundamental Mathematical Equations' ') 


The current equations: 

These equations are sometimes referred to as charge transport 
equations, as they are derived from the Boltzman transport 
equation under the assumption of valid Boltzman statistics. The 


equations relate the electron and hole conduction currents to 
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their drift and diffusion parts and are given as: 


Jn(Xst) = ug(xst) E(xst) nxt + up(x.t)kT BOGE (1) 


(x,t)KT BOGE) (2) 


ca 
— 
x 
ct 
Naat 
! 


p a uy (x,t) ELxet) P(x,t)q = Hy 


J.) = J (xst) - Jn (xt) (3) 


The current continuity equations: 


The equations are statements of conservation of electric charge 


and for this case are written as follows: 


“ti 6a) 
dXeb) af ieee 
aie 1} 2 Jatxet) : 
eet = ot x,t) (5) 


where G(x,t) is the net rate of generation of hole-electron 
pairs due to avalanche multiplication. (recombination has been 
assumed to be negligible) 
The total current equation: 

The total current is expressed as the sum of the electron and 
hole carrier currents plus the displacement current and is 


given as, 


(og Gor) ge = 
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Poisson's equation: 
This equation relates the divergence of the electric field to 
the spacial charge density and can be written in one-dimensional 


form as 


SR Sake 
2 EUG = Bp(x,t) = n(xst) + Nyx) - Nyx) (7) 
Supplementary equations: 


The carrier generation rate!) is defined as, 
G(x,t) = a,(E) v,(xst) n(x.t) + an(E) vp(x,t) p(x.t) (8) 
where the ionization rates a have the form 


Oy p(E) = A exp ( - ( zm ) A, b, m are constants 


These equations are highly nonlinear and to arrive at an analytical 
solution, several simplifying assumptions have to be made, which may limit 
the range of applicability of the resulting solution. Even when using a 
digital computer to solve the equations numerically,it is advantageous to 
make many of these approximations, because of limitations on computer 


memory and to avoid excessively large amounts of computation time. 


Since this study is concerned with transient response, which 
implies solving the equations for a time interval comprising a large number 
of R.F. cycles, to use a minimum number of approximations may be im- 


practical. Therefore, a balance was reached between the amount of simplific- 
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ation of the equations and the cost expended, to obtain a solution 


sufficiently versatile to produce significant results. 


3.3.2 Development of Conduction Current Equation 


The basic equations given earlier, can be simplified with 

very little loss in accuracy by making the following assumptions: 

- the drift current is many times greater than the diffusion 
current so that the diffusion current can be neglected. 

- the electric field is sufficiently high that the carrier 
velocities are approximately constant, at their scattering 
limited values. 

- carrier generation is confined to the avalanche region. 

With these assumptions, the basic equations for the avalanche region 


become , 


J,(x.t) = v, q n(x,t) J (x,t) Shy p(x,t) (1) 
J,(t) = d,(et) + dp(xst) (2) 
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ea) 


acne pete my eat 
notautttb edt nest safest eomtt \ pes 
.bstosfesn 9d neo dnsrwo- roti add ee 02 eas. 
yeiayso std todd dord Uisastoivue at fait statins att 
pntyetisoe halt te , Instenoo Ufedentxoraas 916, eattoofay 


; 


=m 


eouley bottmef 
-Nofpay Sronsisve sd ot banitnos ef notdsisnsp iets. - ‘ 
nofpoy sianslsve oat Tot enotisups stead add .enottqmess — a 
vemoned 
| : 

cr) (tex)q P gv= (tex) ,6 (t.x)o p yy = (tende 
"a 
(S) (Sex) gb +. (Sexdgb-™, (2) gb 9 Wee 
(t ) t ; hn | oat ; * 
X)G 6 . 
+g ae 3 - Gedo ae Ty 


(ty. 


27 


Upon multiplying Eqs. (3) and (4) by q Vi and q v. respectively, the 


p 
result is 
Py J (%st) 9 J,(Xst) 
1 “avn! ape i MpOHtE) Jy(xt) + vyo,(E) d.(x.t) (5) 
en (et) oa (X56) 
Beer Oe eared dla Oh eepabreviiabameeray. 6 


Then,adding Eqs. (5) and (6) and using Fq. (2), the following result is 


obtained: 
Bre) Bee Ret) Dares ts) 
ae ~ Vy oe i Me eee re AN )[o,(E)9, (xt) 3 
GE) (xXst) | (7) 


A more meaningful expression can be obtained by adding and subtracting 


the quantity (v+v )a, (E) J (x,t) in Eq. (7) to give 


p 
> J_(t) eae SAM eye 
a “een “i at a PACT ar eae eles Jo(t) 4 
(a,,(E) - a (E)) Jo(xst)] (8) 


Alternately, by adding and subtracting (v,+V,) a (E) J,(x>t) in Eq. (7) 


the result is 
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ou aCt) 0 J (x,t) eile abd 
mr Pasi a mae Se ee (vitv,) Co (E) J.(t) : 
(a,(E) - an(E)) Jn(x5t)] (9) 


The spatial dependence of Eqs. (8) and (9) can be removed by 


integrating these equations across the avalanche region, with the following 


results; 
a J (t) a Xa 
ae a it J, (Xt) 7 Vy Jpixst J + Crear) elt) fog ED + 
a : 
(vtVp) f (a9(E) = a,(E)) Jy(xst) dx (10) 
O 
2 90(t) Xs ra 
Ky SE = [iy (xt) - Vy Jp(xst)] + (v,t¥p) st) { (F\dx + 
Xa 
(v4v ) f (a, (€) - ap(E)) Jn(xst)dx (11) 


Applying Eq. (10) to the pin junction and using the boundary conditions 
for the carrier currents given in Fig. 3.1, the first term on the right 


hand side of Eq. (10) becomes 
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FIG. 3.1 BOUNDARY CONDITIONS FOR CARRIER CURRENTS IN THE 
AVALANCHE REGION OF (a) ptn JUNCTIONS 
(b) n*tp JUNCTIONS 
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where the relation Je = ne + Jos has been used. 
Similarly, when Eq. (11) is applied to the n’p junction, whose boundary 
conditions are also given in Fig. 3.1, the first term on the right hand 


Side of Eq. (11) becomes the same as that expressed in Ea. (12). 


Using the above results, and after some rearranging of terms, 


the equations become 


Xx 
a 
a J _(t) Vary Vv +v 
S21 (tay ate Cf ato ae n+ RB) as 
— 
chin J_(x,t)d (13) 
( x, ) (a,,(E) - a, (E)) p(x st) dx 
hal am 73 a 
fe) Va, Ve 
Tie ee ) acct) Cf ale) ax 1) (BB) a, + 
O 
Kane 
ee) (14) 
(PY J (ag(E) - a(£)) Jy (est) dx 


where Eq. (13) applies to the pn junctions and Eq.(14) to the n’p junctions. 


These equations can be greatly simplified if the assumption is 
made that a (E) = a(E). This results in both equations being identical 
as the last term on the right hand side of the equations vanishes. Once 
this assumption is made, all distinction between different types of junctions 
and semiconductor materials is lost. Since it has been shown that there 
is a significant difference between the behaviour of n’p and pn brpetione 


and that in general, the ionization rates for holes and electrons are not 
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equal » this assumption will not be made and the last term on the right 


hand side of Eqs. (13) and (14) will be retained. 
3.3.3 Constant Electric Field Model 


In their present form, Eqs. (13) and (14) of the previous sub- 
section are not readily solvable, due to the presence of terms involving 
the integration with respect to the spatial variable (x). One method of 
eliminating these spatial integration terms is to assume a constant 
electric field all across the avalanche region. Then, the ionization rates 


are constant and the spatial dependence of the current densities J,(Xst) 


and Jn (xt) can be taken from the d.c. analysis given by Mo11 16), 
From Moll, the current distributions are given as 
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J J 
a = epee 2 (4) 
ave dee 


for the n’p junction. The term J is actually the instantaneous value 


dic. 
of the conduction current density, J(t), for the constant electric field 
assumption, and using this assumption the current distribution equations 


reduce to 


J (x) e163) 
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for a pn junction and to 
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for a n’p junction. Substituting Eqs. (5) and (6) into Eqs. (2) and (4) 


respectively, these equations become 
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When Eqs. (7) and (8) are substituted into Eqs. (13) and (14) 
(subsection 3.3.2) respectively, and the integration with respect to the 
spatial variable carried out, the differential equations which describe 


the conduction current in the avalanche region become: 


uate )texpl (a, (E)-0 (E))x,J-1} 
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Xa 
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a a, (E e a1, (E) 
(10) 


for pn and n’p junctions, respectively. 


The coefficient of the J(t) term has to obey the criterion for 


breakdown in avalanche diodes, which is given as (16) 
W W 
fence exp[ - i} (a, (E) - a (E)) axa | dx. =e | Gb) 
0 x 


for p'n junctions and (W is width of depletion layer) 
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or 
Xx. = | ve ( o ) (13) 
a a (E)-a_(E n a (E) 


where xX. is the equivalent width of the avalanche region for the constant 
field model. When the electric field across the avalanche region is 

of such magnitude and polarity that the semiconductor junction is just 
beginning to experience significant carrier multiplication, (i.e. the 
breakdown condition has been reached), the coefficient of the J (t) term 
in Eqs. (9) and (10) is zero. Setting this coefficient equal to zero 
results in the same relationship as given in Eq. (13) and Eqs. (9) and 


(10) are consistent with the criterion for breakdown in avalanche diodes. 


3.3.4 Simplified conduction current distribution 


Eqs. (9) and (10), of the previous sub-section, contain a 
large number of exponential terms, as the ionization rates are themselves 
exponential expressions. This makes normalization very difficult and can 
lead to instabilities in the numerical calculations. It would be advantagous 
if alternate expressions for Eqs. (7) and (8), sub-section 3.3.3, could 
be found, which would eliminate the exponential terms in the current 


density equations. 


Alternate expressions for the carrier current densities can be 
found, which are free of exponential terms, by considering the d.c. current 
distribution (16), , the semiconductor junction. These distributions are 


shown in Fig. 3.2 along with an alternate representation that has essentially 
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FIG. 3.2 DISTRIBUTION OF CONDUCTION CURRENTS IN AN 
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the same behaviour. Since the "pte part of the junction has a doping 


19 


level of about 10 cm”? while the "n" part of the junction has a doping 


16 


level of about 10 cm?, the magnitude of Jn (xa) is much greater than 


that of J, (0). (Fig. 3.2). Thus, the relation for the reverse saturation 


current density (J.) can be written as, 


leg diy oa ARG oh es) aE a (1) 


and the alternate representation shown in Fig. 3.2 is valid. 


A further approximation, which leads to very little error, is 


to let J. be zero in the alternate representation shown in Fig. 3.2. The 
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Germanium, which is much less than the typical magnitude of 102cm~2 for 
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magnitude of J, is of the order of 10° "cm 
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n'%a Further, as will be shown in the discussion of the experiments, 


when the Impatt diode is operated with repetitive-pulse-type bias modulation, 
ad.c. current of 1 ma is required to ensure stable operation. This 


reduces the importance of the reverse saturation current. 


The d.c. current distribution for the pn and n’p junctions, 
under the constant field assumption, can be calculated from Eqs. (5), (6), 
(7) and (8) of sub-section 3.3.3. Using the previously discussed 
assumption of Je equal zero, these current distributions are as shown in 


Fig. 3.3. By employing the boundary conditions 


J(0) = J(t) sau 


D (x,/ 8) = J.(t)/2 3d (xa) = 0.0 (2) 
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the current distribution for the pn junction can be empirically expressed 


in the convenient form: 
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J_ (x) 
I(t =3 le 5 ( = i where s>=2 (3) 


Similarly, the current distribution for the np junction can be 


empirically expressed as 


3 Tl (cx de 

ROM ae Care? Ui ee (4) 
Substituting Eqs. (3) and (4) into the conduction current 

relations, Eqs. (13) and (14) sub-section 3.3.2, and carrying out the 


spatial integration, yields the following results: 
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for the pn junction and 
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st here) iii) ve Moe et temp dbsties (BN) lores pc 
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Xa BS) 


for the np junction. 


These latest equations are much simpler than the more exact 
versions given earlier, Eqs. (9) and (10) sub-section 3.3.3, particularly so 
when the terms x, a (E) and x, ay (E) are expanded in a Taylor series as 
powers of the electric field (E). 
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3.3.5 Validity of the Model 


An appreciation of how valid the approximations leading to 
Eqs. (5) and (6), sub-section 3.3.4, are, can be gained by comparing the 
Static behaviour of these equations with the more detailed analysis of 


Schroeder and Haddad! 22) 


nihey eesOlVed GS mull). 9(c)c. (a) pula) sad a eana 
(12) sub-section 3.3.3, to give the spatial dependence of Jp (x) and 

J, (x) for both pn and n’p abrupt junctions. In their analysis of the 
abrupt junction Impatt diode, they assumed a triangular representation for 
the electric field across the junction. This is not entirely accurate, but 
should give just as valid results as the constant field assumption used 

to arrive at the equations derived here. The more exact enaiyet ee )ad vee 
an avalanche region width that is about twice as large for the pn junction 


as for the n’p junction. 


Under the constant field assumption, the relation governing the 
equivalent width of the avalanche region, Eq. (13) sub-section 3.3.3, is 
the same for both pn and n’p junctions. However, if Eqs. (5) and (6) sub- 


section 3.3.3, are solved for the position in the avalanche region where 


n nt i p aver 
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m ] 
> | Seep) ) #9 Ca cepa ey (2) 


for np junctions. With reference to Fig. 3.3 and Eq. (13) sub-section 


3.3.3, the following ratios hold, 


xy sees tb a (E) +0 (E) a (E) 

Ke = i gn ( 5 aE) ) in ( ay (E) ) (3) 
for pn junctions and 

x, fe 2a,(E) a, (E) 

a s .* an ( a, (E)¥a, (6) y/ an 4 ay E ) (4) 


for n'p junctions. These last relationships can be evaluated by using 


(14) 


data given by Sze for the ionization rates of Silicon where 


a (E) = 3.8x10° exp(-1.75x10°/E) and a (E)= 2.25x10/ exp(-3.26x10°/E) 
(5) 


Using a similar procedure to that outlined in Eqs. (1) to (4), 
values of X374/Xq and X9710/%q may be calculated. This information, along 
with other pertinant data, is summarized in Table 3.1, for the constant 
field model. The point at which the current levels reach 90% of their 
final values may be taken as the extent of the equivalent avalanche region 
for the model. Table 3.1 shows that the use of the constant field 
assumption gives an avalanche region width for p'n junctions that is 1.54 


to 1.3 times as large as that for the np junctions. This is in general 
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TABLE 3.1 


SUMMARY OF PARAMETERS FOR CONSTANT FIELD MODEL (SILICON) 


PARAMETER ELECIRGCe RLELD ( x10° V/Cm ) 
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agreement with the previous work 22), 


Previously published Wonkeeee have produced many graphs 
showing the predicted values of the avalanche region width, depletion 
region width and, breakdown voltage with reference to the doping levels of the 
semiconductor. Unfortunately, the actual diodes do not closely obey these 
theoretical predictions. The actual diodes generally have a depletion 
region width that is less than the predicted value, while the avalanche 
region width and breakdown voltage are greater than the predicted value. 
This restricts the applicability of detailed calculations when analysing 


Impatt diodes. 


The discrepancy between the actual diodes and the detailed 
calculations is due, in part, to the difference between the actual and the 
theoretical distribution of the electric field across the diode junctions. 
This can be observed by comparing the electric field,used by erhamrete en ee: 
which closely resembles that in an actual diode, with those used in the 
various theories. Fig. 3.4 shows a reproduction of the doping profile and 


(24) 


electric field distribution used by Scharfetter The value of Xs given 
in Fig. 3.4 was calculated by applying the given electric field to the 
ionization integral of Eq. (11), sub-section 3.3.3, where the equation was 
evaluated to the 95% point. Also included in this Figure are the equivalent 
constant field and triangular representations for the electric field. 

These equivalent field distributions give the same voltage drops across the 
avalanche and drift regions as does the given electric field. It would 


appear from Fig. 3.4 that the constant field model is equally as valid as 


the triangular field approximation. 
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From the preceeding arguments, it is seen that the use of the 
constant field model retains a measure of the difference between pn and 
np junctions for Impatt diodes and also presents an equally valid 
alternative to the more commonly used triangular field representation, as 
shown in Fig. 3.4. Further, since detailed calculations involving doping 
densities etc., have some shortcomings when describing the behaviour of 
actual diodes, the flexibility inherent in the representation of the constant 


field model, shown in Fig. 3.3, is very attractive. 


3.3.6 Development of Electric Field Equation for Avalanche Region 


So far, an expression for the conduction current in the 
avalanche region has been formulated as a function of the electric field. 
It still remains to derive an expression for the electric field across the 


avalanche region. This required equation can be obtained by considering 


Poisson's equation (25) 
3 Ext) = a (Ny - Ny + p(xst) - n(x,t)) (1) 
S 


which, when rearranged and integrated becomes 
X X 
E(x stpe fc Ny SuNgiichxttse ch fore - n(x,t))dx + E(o,t) (2) 
5G Ss 
fe) fe) 


or 


E(x,t) = E-(x) + SPACE CHARGE + E(o,t) 


where the zero value of the spatial variable (x) is as shown in Fig. 3.4 
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The term Ee (x) is the electric field due to the fixed charges and is 
independent of time. The space charge term is the electric field due to 
the conduction current, while the term E(o,t) is an arbitrary constant for 
each interval of time. This "constant" term consists of a d.c. and an 

a.c. part: the d.c. part from the bias voltage applied above breakdown 

and the a.c. part from the R.F. voltage applied across the diode. 

This is more readily apparent when the Eq. (2) is rearranged and integrated 


to give the voltage across the diode. 


Rewriting Eq. (2) as 
Xx 
E(o,t) = E(x,t) - E(x) - a J (eGot)-n(xt) a (3) 
° 6) 


and integrating across the diode, whose depletion width is W, results in 


a voltage relationship: 


W W X 
E(ost) W= V(xst) = fepx) - f Lf (rOct)-n(xst))dxt ax (4) 
fe) 6) 


226 
where the term, V(x,t), is the total voltage across the diode and consists 
of the voltage drop due to the fixed charges of the junction, the additional 
voltage to bias the diode above breakdown, plus the R.F. voltage across 
the diode junction. Taking the algebraic sum of the first two terms on the 
right hand side of Eq. (4) yields a d.c. and an a.c. term, with the d.c. term 
being essentially the bias voltage above breakdown. Keeping the above 


arguments in mind, Eq. (4) is rewritten as 
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W xX 
E(o,t) = rT? 7 i fe j (p(x,t) - n(x,t))dx'dx (5) 
26 0 


where Vp is the bias voltage above breakdown and ve Z is the R.F. volt- 


age across the semiconductor junction. 


Eq. (5) is of the same form as that initially derived by 
eae used in practically all discussions involving Impatt diodes. The 
last term on the right hand side is due to the space charge of the 
conduction current and is usually only applied to the drift region of the 
diode, as it is assumed that the avalanche region is very narrow. The 
assumption of an avalanche region of negligable width is not valid for most 


(26) that the avalanche region can occupy 


actual diodes as it has been shown, 
an appreciable portion of the total depletion region of the diode. Also, 
the voltage drop across the avalanche region is approximately ' the 


(24) Consequently, there should be a 


total voltage drop across the diode 
significant space charge effect from the conduction current in the 


avalanche region. 


With the constant field assumption and using the d.c. current 
distributions given earlier, it is possible to apply the space charge 
term, independently, to both the avalanche region and the drift region. 
Rewriting the space charge term of Eq. (5) as a function of currents, it 


becomes, for the avalanche region 
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Making use of the current density relations given earlier in sub-section 


3.3.4, Eq. (3) and repeated here for reference purposes, 


uate) 1 b Jo(x st) J_(x,t) 
ey one Jrsgs * ei x, ) 78) AL I(t) 


where, these relations apply to a p'n junction. Substituting the above 


relations into Eq. (6), the result is 


ey ae Vn he Aaa 
SPAGE GHARGE (GAVAle) =. | —4 -f J (€) ee) 
W € C VV 2 V 
ae np n 
( $X yP 4 ax! ax (7) 
xX 
a 
and after carrying out the integration 
V xe 
all Dye bead aL “1 
SPACE CHARGE (AVAL) = We. inp Lee ) ee ] d(t) 


(8) 
which is valid for a p'n junction. For an n’p junction simply replace "b" 


by "d" as defined previously in sub-section 3.3.4. £q. (4). 


In applying the space charge term to the drift region, use is 
made of the assumption that avalanche multiplication does not take place 
in the drift region. In the drift region, the conduction current consists 
of bunched minority carriers. For a single type of carrier in the non 


. : +, , 
dispersive drift region the following relations are valid for a p n junction, 


svods art a a smofsonut a ty 5 oF Y i in a 


at tfueey " (3) lata a 


(ta gS ae | Ny noc 


? ; i) 
(s) xb 'xb [ “( a aii 5 
ste 
notteypstnt oft do eniyyTEs: sods bas 
qn’ I at = (4AVA) 3 
(3)6  copayrteay | ~V Jag AG oe ( JAVA) nia At 
(8) . Ts 


"d" 9961997 xiqate noftonut q° Ans vod .norsonul n’q 5 ToT bhTev 2 
(Bh) .pa .B.£.€ nottose-due nf yfevotverg banttsb os hs 


2h seu ,notpe ttinb ons ot mist sprsrio sosqe2 ery ontyfags. al 
snag ast Jon “ot notseatfqts Tum sionsieve Sets nottqmess: ort to 9 
exetenoo, nevus fot toubno> oft .notpet ttith ond mi a SI ng 
aon oa ni i9ty69 to sqyt afonte B 10% .219hr169 yd tronfm barton 
enobtonut 9 +g 6 10? bilev ot 2nottefor antwor ror ant rrr yy ta 


, 
—< 


48 


The space charge term of Eq. (5), when applied to the drift region 


becomes 


WwW. W 
SPACE CHARGE (DRIFT) = fz ice (x,t))dx!' dx (9) 
S 
Xx 


aa 
where only one type of carrier is permitted. Continuing the derivation 


for a pn junction and carrying out the first integration gives; 


SPACE CHARGE (DRIFT) = Wes {ite x) ae Ce a de (10) 
(@) n n 


where W'is the width of the drift region and the relationship for a single 
type of carrier is used. This last equation can be further simplified by 


a change of variables, such as letting 
ber 2 
RAM ae 


Using this change of variables the space charge term becomes 


< 


t 
SPACE CHARGE (DRIFT) = We | (Ty-t+t') J.(t')dt' (11) 
S 
t-Ty 
where Ty is the transit time for the carriers in the drift region. 
Substituting Eqs. (8) and (11) into Eq. (5), the relationship 


for the electric field in the avalanche region becomes; 
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which is valid for pn junctions when « is set to unity. (For n’p 


junctions substitute "d" for "b" and « is set to beg Sel 


It is interesting to note that the term relating to the space 
charge in the avalanche region varies as the square of the equivalent 
width of the avalanche region. Otherwise, the space charge term for the 
avalanche region behaves as expected, anere its magnitude is directly 
proportional to the magnitude of the conduction current and vanishes for 


small values of Xo: 


For n’p junctions, the space charge associated with the avalanche 
region will be significantly less, approximately a factor of 4, as the 
value of Xs is approximately % that of pin junctions. In the constant 
field model, the equivalent avalanche region widths for np junctions are 
approximately 2/3 those of the p'n junctions. (reference Table 3.1) 
However, the magnitude of the "d" term for np junctions is less than the 
magnitude of the "b" term for pn junctions. This further reduces the 
value of the avalanche region space charge term for n’p junctions and Eq. 
(12) should adequately describe the difference between the two different 


types of semiconductor diode junctions. 


3.4 Summary 

The foregoing analysis has been carried out with the pre- 
reauisite that the equations developed would be used for transient analysis 
and would be solved on a computer. Therefore, considering the large 
number of R.F. cycles necessary to characterize transient response, the 
equations should be sufficiently simple that the solutions not take an 


excessive amount of computer time. The following equations that describe 
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the conduction current and the electric field in the Impatt diode are 
the results of the foregoing analysis and are summarised here for 


reference purposes. 


The conduction current equation for the pn, abrupt junction, 


Impatt diode is 


1 
Reet se) 


while the equation describing the electric field in this semiconductor 


junction is 


2 
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The equations describing the np junction are identical in form but with 


appropriate changes in some coefficients. 


These equations satisfy the preceeding criterion of simplicity 
and ease of solution, while at the same time using a minimum of approximat- 
: + 
ions. The equations are sufficiently versatile to describe either p n or 


+, ; : : . ? 
np junctions in a wide variety of semiconductor materials. The accuracy 
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of the equations should be adequate, as the assumptions used in 
developing the current equations have been shown to be valid and the 
electric field equation includes an additional term to account for the 


space charge effect of the avalanche region. 
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CHAPTER IV 
COMPUTER SOLUTIONS OF THE ANALYTICAL EQUATIONS 


4.1 Introduction 


In this chapter, the analytical equations derived previously 
in Chapter III, are solved. Two approaches are discussed, as both Analog 
and Digital computers were used. Computer usage provides a fast and 
accurate method of obtaining the equation solutions. Thus, a study of 


the results for a variety of parameter values can be readily carried out. 


Obtaining solutions for a range of parameters is very important 
when comparing the computations to experimental results. There are 
parameters of the Impatt diode that cannot be determined without some 
trial-and-error. Although the depletion region width and doping con- 
centrations can be estimated from capacitance - voltage measurements, 
the division of this depleted region into separate avalanche and drift 


regions is somewhat arbitrary. 


Since many computer runs needed to be carried out, it was thought 
to be advantageous to normalize the analytical equations so that the 
required accuracy could be attained using the simplest computational tech- 
niques. For digital computers, this would be single-precision arithmetic, 
whose usage would reduce the computation time and computer memory space 


required for each program run. 


4.2 Normalization of the Analytical Equations 
Equations describing the Impatt diode that were derived in 
Chapter III are repeated here for reference purposes. The current equation 
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ad (t) vtv 1 b 
xt = ( 3 J xt Jt) xX. [ a, (E) ( mye cf ay (E) ( ay? 3 
14a, " 
a 
and the electric field across the avalanche region is 
Vp Va Cc Vy x, ] Yn “p ] 
Bloat * ay EW yaa not as V Tamra ce thle 
t 
’ f (I, t+ t') a(t!) dt) (2) 
t-Ty 


where « is unity for the pn junction. {ror np junctions, the equations 
are the same except "b" is replaced by "d" and a, (E) and a ,(E) ineedem 1) 
are interchanged. Also, in Eq. (2), x = Vo Mn for the n’p junction. 

To facilitate the use of compatible magnitudes in the computer 
solutions, the above equations were time-scaled and rewritten in terms of 


currents. The time normalization was carried out according to the relation, 
c= 10% (3) 


where t is the scaled time and t is the actual time. The equations are now 


rewritten in terms of t and using the relation 


Ect )¥= foo dA = J_(t) A (4) 


where A is the effective cross-sectional area of the p-n junction. The 


time normalized equations in terms of currents are 
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dd _(t) Vv +V 
Cc -10 n 
een ee ee ry ee ar 
- ] + ie (5) 
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Using the values for ve and Vy that are given in the references 


where 


v=10/ cm/sec and Vo Sep > & 107 cm/sec 


and also rescaling other quantities as follows: 
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Before the current equation is further normalized, it is 
advantageous to expand the ionization rates in a Taylor series about the 
variable (ESE e)s where E. represents the value of the electric field at 
breakdown. The resulting relations for a (E) and a (E) will be identical 
except for notational differences, so only the a, (E) equation will be 
fully developed. 


The ionization rates are given py (14) 


a, (E) = F exp (- f/E) and a, (E) = G exp ( - g/E) 


where 
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Expanding a, (E) in a Taylor series and taking only the first three terms, 
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Eq. (8)is a rapidly decreasing function of MES and the third term 
of the variable contributes very little to the calculation, but will 


be carried for completeness. Similarly, the expression for a, (E) is 


-baqofsvab wt 7 
(Dy novig avs eater notaastaoh aa 
(3\p - ) axa B= (B)go bas (Bt -) axe 


4 Ole 
aia 

‘ore assea ; MOreeeeY | 
Sor x ase *p  : Oh Nt =F | 


ames sondt gett sft vino pntdet brs e9trae yofyst 5 nt (3), ef 


(a) nee S(_a=a) (2) 0 rant 
ic a tgs en (Ae er 
) 2 | E> 4B, aa wag 
oe il ; See vi - 


5373 O36 lh 
te tate Gia tt eee “f 
(8) (i 3-3) ie 2° Be A oe or 


mist byt? sdt brs (.3-3) to alae ann Po 
Titw Jud ,notsetustss odd of 9fttl yey zedudtysnos ofdstrsy end 
2h (3) go 10% nohezsygxe ott list : f __ raeanenatignae 


56 


When Eqs. (8) and (9) are added, with the appropriate coefficients 


applied, the inner bracketed term of Eq. (5) becomes, 


x, (By) aplEg) - 11 +0 x, ( pat) anle) a4 
¢ 
x, (Py) a (EQ) se! (E-E.) +... (10) 
Cc 


The first term on the right hand side of the above equation is the break- 
down condition discussed in sub-section 3.3.3. Therefore, this term 


vanishes from the equation and Eq. (5) can be written as 
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With reference to Eq. (5) of sub-section 3.3.6, the variable 
poe) is the same as E(o,t) defined in Eq. (6) and (7) here. By applying 


the normalization factors used in Eq. (7), Eq. (11) becomes 


respectively, 
nora = °2. (E) frthe + [ (i é fe th youn 
Jenlai alas Split eaiaan OT 
E 5 2 
C c Cc 
[ x ( = Ne - t+ — E(o,t)” } (15) 
E. EY EY 
and 
mise ny b Pulse. genes 1 eto, ovnare tee 
BETA = 2 ( erp a aE.) { ea ECO.) rel x I soa = salut) "3 
C Cc c 
Nn) aa g 
[Le(4y)> - 46+ 4q] Elo.t)” } (16) 
oo ie go ee 


4 
4 : nowy a . 

; 
: 7 - a 

a u 


iy : 


ee) ed + CDS Nee 
ee es 


nh 


+ (3-3) [¢= - 
32 
(ef) (3-3) Sree “gy 


- 


5. 


i ai 
aldsfrév aft .8.€,€ norsoase-due to (2) .pa of sonieysTs" agtw 
onrylags yo °.ayert ({) bas (8) .pa at bsniteb f20)3 26 ams oat 2 a 


zamooed (Tf) .p3 .(X) .pa nf beew arot58 notssst lennon 8 (A: 


. 
(7) Ts 
(AF) ( T+ ( ARG. + WATIA] (7) 51) e.Ni = ae 


pnttise vou enw 


“OE (3 go lA ge «Bete ala ( ge xf aye . tag 


(ef) { 


i 


58 


The magnitudes of E. and x, are related through the constant field approx- 
imation where EO is the value of the electric field in the avalanche 

region, for a given value of the avalanche region width, Xa This relation- 
Ship is formulated in sub-section 3.3.3, Eq. (13) and also summarized in 
Table 3.1. As an indication of the magnitudes that can be expected for 

the coefficients of Eq. (14), when the following representative values of 


the variables are used 


mI 
" 


Sonos a (E.) Se O7 Sy a (E.) =n). Wy Xie 173 


pe a2. 10st 1 59g = 26 


Eq. (14) for the current in the avalanche region becomes 


——— = oa T(t) [0.173 Efo,r) + 0.012 Efo,x)* + 0.003 Efo,x) HT, } 
a 
Eqs. (14), (15), (16) and (7) describe the conduction current in 
the avalanche region of the Impatt diode in a form that is very readily 
computed. However, further equations are necessary to relate this avalanche 
current to the current flow in the circuit external to the semiconductor 
diode. This external conduction current a can be formulated by 
considering that a pulse of carriers is formed in the avalanche region and 
then "drifts" across the drift region of the diode without dispersing. 
As this current pulse travels across the drift region, whose transit time 
1ShT 
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Causes a current to flow in the external circuit. This external conduction 


current can be written as 
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after magnitude and time scaling. 


The total current (I~) in the external circuit of the semiconductor 


-) 
diode is the sum of the conduction and displacement currents, which can 


be written, after magnitude and time scaling,as 


Eis teem Ao (18) 


d OT 


where Ca is the depletion region capacitance in picofarads and V(t) is 


the R.F. Voltage across the diode junction. 


In the computer program, the Impatt diode is completely described 
by Eqs. (14), (17), (18) and a relation for V(t). V(t) can be expressed 
as a simple trigonometric relation such as a cosine or sine function. The 
waveforms given by Eq. (18) and the relation for V(t) are Fourier-analyzed 
and the complex ratio of the fundamental frequency components taken to 
obtain the electronic impedance We) of the semiconductor diode. The 


Fourier analysis procedure is explained in Appendix A. 


The effects of temperature changes upon the diode junction are 
discussed in the following Chapter, where it is shown that this effect 
can be taken into account by adjusting the value of Vp insEGua/JeaSed 


function of time. It is also shown in Chapter V, that, for the short- 
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duration bias voltage pulses used here for the transient measurements, 
the effects of junction heating are minimal. Therefore, the effects of 


the junction temperature changes will be ignored in this chapter. 


4.3 Analog Computation 


Initially, the normalized equations were programmed on the 
Pace" analogue Computer. This served a two-fold purpose; (1) to gain 
familiarity with the behavior of the model and (2) the possibility of 
constructing an analog model to describe the transient behavior appeared 
promising, as steady-state models of the Impatt diode have been programmed 


on analog computers (28) 


Analog computers have a limited number of amplifiers and there- 
fore, the extent of the non-linearity that can be programmed is restricted. 
Hence, any representation of the Impatt diode that was modelled on the 
analog computer would need to be simplified. Even so, an approximate analog 
model of the Impatt diode would be valuable for demonstrative purposes, so 


work was carried out towards this goal. 


4.3.1 Equations for the Analog Model 


The model of the Impatt diode used for analog computation was 
essentially the same as that given by the equations derived in Section 4.2, 
Eqs. (14), (7) and (18). The simplification was mainly made in the 
current equation, where only the first term in E(o,t) was taken. Also, the 


model was assumed to have equal ionization rates and drift velocities for 


* 
PACE - Precision Analog Computer Equipment; a general purpose d.c. 


analog computer, manufactured by Electronic Associates Limited 


(Model 231R). 
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both electrons and holes. The resulting equation for the current in the 


avalanche region was thus, 


ee GREY tay te ie (1) 
a a 


where Ta is the transit time of the avalanche region. 


The equation for the electric field is the same as that given 
by Eq. (7) in Section 4.2, except that the space-charge term associated with 
the avalanche region was neglected. To facilitate programming of this 


equation, it was written in differential form as 


3 Elo, foal oe atl 1 SCH ASi 
ee he = pf tylgt) - (Te )de'] (2) 
Sad 
Ly, 


where the bracketed term originates from the relation 


T ih 
2 fern) Te(re ar were Ure) - | ates )dr* (3) 
ar ToT 


and a d.c. initial condition simulates the constants in the original 
equation. Similarly, the expression for the external conduction current 


becomes 


] bs Je 
ts t [ I(t) - Tei) ] (4) 


The R.F. voltage across the semiconductor diode was simulated 


by setting 


i 
. 
i 


» ¢é ; fy 

a 
notps? sdonslsvs sft to omit cudnt iikiog 

s ; — : 7 a 

navie. tedt 25 omse oft et bioft otvdoels odd vot notdsups oT i ah 
dgtw betetooez6 aie? Spvedo-sonqe odt fant tqvoxe ,5.* notsose at (X) pad 
athit to pntimsipoyg stair {fost of -batosipsn 2sw ao tps orionsteve oA : 
2 mot fettnovsitth at netttw zew tt “not taupe - | 


75 ’ ; 
3) | prsbt' 2) (i tek es ae 
b2 W . 
tt \» | 


> 
nofislay eds mov? estantotto anes betstoerd wae 


z rt ; 
(t) 'sb( x) a1 \ (s)GT ype tb CGT Cres ~f - 
b' « a i 


fentetvo ott nt 2tnsi2znoa sit estslumte nots rbnos isttint .o.b 6 brs a 
Insviva noksaubnoo (envedxs ahd ot nofezsiqxs ads ,yVltsrtmi2 «mort URo =] 


- aomozed 


(d) | [ (arent “ () 7 - 
beteiumie ew sboth voioubnostmez 9d 220798 apatfoy .4.s sit. 


1g 


_" Fa nvelieay 


62 


V(t) = M Cos (w rt) (5) 
where M is the amplitude constant; the programmable differential equation 


for the R.F. voltage then becomes 


2 
Che Aa Ver) (6) 


Since the displacement current is directly proportional to the time 
derivative of the R.F. voltage, it was analog modelled by sampling an 
appropriate portion of the first derivative signal in the analog program 


for Eq. \io). 


A schematic diagram of the program used is shown in Fig. 4.1. 
This program contains some extra components, such as absolute value circuits, 
comparators and redundant amplifiers, all of which were found necessary 
to improve the stability of the program. This is a fairly complex program 
and practically al} the available amplifiers and integrators of this 


quite large analog computer were used. 
4.3.2 Results of the Analog Computations 


Steady-state operation of the model was successfully obtained. 
However, some difficulties were encountered in obtaining solutions to the 
equations and the model was found to be too unstable to produce accurate 
results under transient conditions. Some of the instabilities were inherent 
in the programming, but the main causes of the unstable results were 


associated with equipment limitation. 


One of the inherent weaknesses of the program lay in the sim- 


ulation of the time delay necessary to model Eq. (4) of sub-section 4.3.1. 
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FIG. 4.1 ANALOG COMPUTER PROGRAM FOR AN IMPATT DIODE OSCILLATOR 
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The circuit used for the time delay is shown in Fig. 4.2. 


(29) 


This is a Pade approximation of the fourth order. Although this 
circuit is among the most accurate approximations that can be made on 
the analog computer, it still suffers from d.c. drift and ringing when 
a sharp current pulse travels through it. One way of overcoming these 
problems would be by simulating time delay, by using a digital computer 
for storage, with analog-to-digital and digital-to-analog converter 


equipment. As the required converter equipment was not available, this 


corrective action was not taken in this study. 


The major equipment limitations were the d.c. drift of the 
integrators as well as noise in the output of the multiplier. Absolute 
value circuits and comparators were placed after the multiplier and 
in the center of the time delay circuit. Although these extra circuits 
reduced these sources of error, they did not eliminate them. Consequently, 
the transient results were not very accurate and were difficult to 


reproduce. 


The above - mentioned problems were not as severe when a strong 
R.F. voltage was applied to the Impatt model and the system operated under 
steady-state conditions. However, the accuracy and reproducability of 


the results were still somewhat limited. 


4.4 Solutions by Digital Computer 


An analysis of the transient behaviour of the Impatt diode re- 
quires a solution over a large number of R.F. cycles; thus, each computer 


run may be of appreciable length. Also, to reconcile computer behaviour 
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of the Impatt diode with the experimental results, some adjustment of 

the diode's parameters will have to be made, and this may necessitate 
running the program several times. It would be advantageous to use 
Single-precision arithmetic, as this allows a reduction in the computation 


time and the amount of computer memory required for each program run. 


By casting the equations, describing the characterization of 
the Impatt diode, into normalized form, it is possible to use single- 


precision arithmetic and still maintain sufficient accuracy. 


The normalized equations derived earlier, Eqs. (7), (14) and 
(18) of Section 4.2, were programmed on an IBM 360/67 digital computer. 
A further advantage of using the normalized version of the describing 
equations is that the resultant solutions are very stable. Consequently, 
less stringent error critera and larger time intervals can be used in 


the solution of the equations, while still maintaining sufficient accuracy. 


4.4.1 Procedure for Numerical Analysis 


The set of equations to be solved include only one differential 
equation; however, if a coupled cavity were included in the model, there 
would be three first-order differential equations that would require 
Simultaneous solutions. To ensure flexibility, a routine was used that 
would allow the solution of one-or-several differential equations. A 
specific subroutine for the solution of the differential equations was 
not written, as several excellent subroutines are already in existence 


and little would be gained by duplicating the existing programs. 


A subroutine employing Hamming's modified predictor-corrector 
method, with a fourth-order Runga Kutta method to supply the starting 
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This subroutine is faster than the fourth-order Runga Kutta method 

alone and has a better error criterion. The program for the solution 

of the set of equations was built around this HPCG subroutine, as it 

had separate input and output subroutines of its own, that were supplied 


by the user. 


In programming the normalized equations for the Impatt diode 

(Eqs. (7), (14), (18) of Section 4.2 plus a relation for the R.F. voltage) 
the equation describing the electric field in the avalanche region was 
rearranged slightly; the form actually used in the program was expressed 
as follows: 

(VitVe) Vo (t) . 

ep ee ral cata ay eT ae 

W W 1.04WA 
AGES 
d 
(1) 
mf 

Where V, is the d.c. bias voltage, above breakdown, needed to maintain 
a given steady-state current present during the experimental measurements. 


The term V, still applies to the bias voltage pulse applied after the 


B 
reverse breakdown voltage has been reached. A further change that is 
evident in Eq. (1) is that the coefficient of the space-charge term per- 
taining to the avalanche region has been condensed into one constant, 


K+ The latter change made this coefficient more readily available for 


variation when comparing the computed and experimental results. 


The relation describing the R.F. voltage across the semiconductor 


diode (V(t)) was chosen as 


V(t) = M Cos (w t) (2) 
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where M is an assigned magnitude and w is normalized in accordance 
with time scaling relationship for t (Eq. (3) Section 4.2). Here, 


the R.F. voltage is assumed to be at a single frequency. 


An important problem encountered in the use of this program, 
was in determining the most practical size of time increment for the 
solution of the differential equation. It was readily determined that 
a time increment of 0.01, for the normalized time t, would produce 
practically the same results as any smaller time increment used. How- 
ever, a problem developed in that the increment size would gradually 
change as the program time progressed. By using the time scaling of 


10/9 


» the program time variable has to run to 300 to simulate 30 nano- 
seconds of actual time. In so doing, the change in the effective size 
of the time increment is sufficient to cause significant inaccuracies 


in the Fourier analysis results. 


The solution to this problem lies in choosing a size of time 
increment that is convertible exactly into hexadecimal so that all the 
Significance is in the first few digits. The computer stores all inform- 
ation in hexadecimal and the program increases its time variable by 
continually adding the increment size to the existing time. As the 
existing time becomes greater than 16, during the process of addition 
some of the significance of the time increment is lost. For values of 
the existing time greater than 256, even more of the significance of the 
time increment is lost during the addition process. Consequently, the 


effective size of the time increment appears to decrease as the total 


empypo1g 2tdt to seu sit nf BSeenUANs violate ry nit a : 
iz 10? tnemevont omft to este Teots551q 320m ant gntatmrsteb at 26H | 
sef9 bentavaysb yfibsey 26W +1 .nofssups ietgnatsttib oft To notdutoz *% oar 
soubovq biuaw ,r amtt bestlemyon sit rot , 10.0 Yo dnamerant smt3 en 
-woh .bseu trevietont omrt yollsme va 26 23 vest omée ant vifestaoeyq 


oe? ae 
uffeubsip bluow sst2 tnemsvont say Jedd nt baqoloveb nat dong 5 Pi ; 


to pnifsoe ombs att onteu y8 .bseesipotg smtt mevporg sid 26 ae : 
-onéy OF Sisfumte of OOE ot muy oF ead sfdetrsy omty mayporg oy Mor 7 
este avtiastte oft nt epredo etd .pntob o2 ni. .smtt Taytss To zbnoee ve 
satosiuoosnt tneatttapte s2ueo of Jnotorttue 2t snpmstont smtt sft to Ve 
.2tfuesy etaylens watwo? addonk | 
ant# to ostz 8 entacoto nt zatlimatdotg edt of noteifor SAT chs * 
ont {fe tent oz Toutssbexed osnt Yftosxo sidttyvsvnoo ef Jedd nemetont 
-miotat [Ts 2a%ot2 votugmoa sAT .2tipib wot sevtt odd nt ef nites 
ud ofdetrev omtd ett egebsiont meveorg ond bys femtosbsxed nt norte 
oat 2A .ombt onttetxe ont ot asta Snemetsnt odd pribbs iF eunt nos 1 
notatobs Yo 2299070 Sad pntiwb .a! nsrit yIso vp zamonsd smtz onisetns 
to vaufev 101 .teol 2t dromatoat smtt sad Yo sonsorttngte Say Yo onioz . 
ant Yo sonsoftinete edt to s10m Nave. ,de8 apdt issesvp ons webiste eee 
ont .y¥ftnsupsenod .229o001q nofttbbs- oni pentru 201 et 


lefot sit 26 se6o1D9b OF iene Seaiaa smtt oz to este witostts 


ay 
‘ we 


69 


scaled time increases. However, if the size of the increment is 
properly chosen, a value of 1/128 was found to function properly, there 
is no loss of significance during the progress of the program time 


and the effective size of the time increment remains constant. 


The selection of such odd sizes of the time increment as 1/128 
may necessitate some adjustment in the number of points taken for the 
Fourier analysis. However, any problems such as this can usually be 
alleviated by a judicious choice of frequency for the R.F. voltage 


across the diode. 


4.4.2 Results of Computations 


In programming the normalized equations, values need to be 
assigned to several parameters. These parameters are listed in Table 4.1 
and the values assigned are those used in the final comparison of the 
computed and experimental results. The depletion width (W) was obtained 
from capacitance-voltage measurements taken at | MHz and the cross- 
sectional area (A) was obtained by examining a segmented diode package 


under a high power microscope. 


During the analysis of the computer results it became evident 
that the equations were more sensitive to some parameters than others. 
The most critical parameters were Tq? C yok and Xa Small variations in 
W or A caused relatively small changes in the computed results. While 
the equations were more sensitive to variations in K and x they were 


most sensitive to variations in Tq and Ca 
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TABLE 4.1 


PARAMETER VALUES FOR THE NORMALIZED EQUATIONS 


T iC 

d d 
1.451 0.26 |0.36 
1.45 | 0.26 |0.36 


WeOSHV0 e327 710530 
i Soler) WWE Syer/ a iia) 


The value of the carrier transit time for the drift region 


(ty) has a strong effect upon the computed results. It is this parameter 
that largely governs the d.c. current response of the equations to a 
given bias voltage. Increasing the value of Tq decreases the level of 
d.c. current in the simulated Impatt diode, for a given value of bias 
voltage. Consequently, an estimate of the value of Ty can be 

obtained from a comparison of the current levels in the simulated and 
actual diodes for a given bias voltage. A better estimate of the 

value of Tq can be obtained by determining the width of the drift 
region and dividing this by the scattering-limited velocity of the 
carrier‘ #4) An estimate of the drift region width can be obtained from 
extrapolating the value of Xs used in the program, back to an actual 
width of the avalanche region and then taking the difference between 
this actual width and the value of W. Once determined, the value of 


Tq is held constant for a particular diode. 


The depletion layer capacitance (Cy) greatly influences the 


magnitude of the calculated electronic impedance for the semiconductor 
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diode junction. This is well illustrated in Figs. 4.3 and 4.4, 

where the only difference is in the value of Cas which was 0.37 for 
Fig. 4.3 and 0.35 for Fig. 4.4. It is interesting to note that the 
influence of this parameter is mainly in its effect upon the magnitude 
of the computed results and not upon the general shape of the plotted 
curve. This parameter was initially set in accordance with the value 
found from the 1 MHz capacitance-voltage measurements and then 
adjusted slightly to obtain a better comparison between computed and 


measured results. 


The main influence of the parameter K is in the slope and 
curvature of the computed plot of the electronic impedance, while 
changes in magnitudes of the plotted values are influenced to a lesser 
degree. This effect is illustrated in Fig. 4.5, where the only 
difference in the parameter values is in the magnitude of K. The 
time or value of d.c. current where the electronic impedance plot 
“changes direction" is a sensitive function of the value of K. This 
parameter is initially estimated from the theory and then adjusted to 


correlate the computed and experimental results. 


A value for the equivalent width of the avalanche region (x,) 
can be arrived at by considering Figs. 4.4 and 4.6. The value selected 
for Xo is seen to have a distinct effect upon the plotted characteristics 
of the electronic impedance. For the higher value of x, (Fig. 4.6) the 
curvature of the impedance plot is reversed in the region of 7 nanosec- 
onds and, also, the magnitude of the plotted values in this region 


are reduced somewhat. From the characterization of Scharfetter's model, 
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Fig. 3.4, the equivalent width of the avalanche region is known for 


om 3. An estimate for the doping 


a doping level of about 1.1 x 10 
levels in the diode junction can be obtained from the 1 MHz capacitance- 
voltage measurements and then the above information can be used in 


conjunction with the experimental results to determine the value of x: 


In Figs. 4.3 to 4.6 the values of the parameters were for 
diode #1, driven by a R.F. voltage whose frequency was 6.0 GHz. The 
values of W and Tq were fixed at 3.75 and 0.26 respectively, with the 
value of K set at 0.0238 unless otherwise specified. The same bias 


pulse was used for Figs. 4.3, 4.4 and 4.6. 


During the computations, it became evident that, for a 
given value of bias voltage pulse, the magnitude of the d.c. current 
response observed in the computer model of the Impatt diode, was not 
quite the same as that observed in the experimental measurements. In 
diode #1 the computed current response had a greater magnitude than 
that experimentally observed for a given value of bias voltage pulse, 
while in diode #8 the opposite effect was found. This discrepancy 
could not be explained by a choice of the parameter values in the 
computer model and was thought to be a result of the manufacturing 
techniques used in fabricating the Impatt diode. Therefore, all the 
computed results were compared with the experiments on the basis of 
their d.c. current response and not on the exact magnitude of the bias 


voltage used. 


The final computed results were plotted along with their 


associated experimental results given in Chapter VI. The series 
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resistance (R.) due to contact losses, the resistance of the substrate 
and the undepleted region of the diode has not been taken into account. 
However, as the computations give the electronic impedance, the value 
of R. can simply be added algebraically to the computed results. 

This action will result in a shifting of the plotted computed results 


along the real axis of the graphs. 


4.4.3 Estimation of Computational Errors 


The error in numerical calculations consists of truncation 
(discretization) and round-off errors. These errors are mainly 
controlled by the type of method used to solve the differential equation(s) 
and the selection of the step size. In choosing the step size, a 
balance must be maintained between the accuracy that is acceptable 


and the time required to complete a program run. 


In general, the subject of round-off-error propagation 
is poorly Carne accumulated round-off-error is not simply 
the sum of the local round-off-errors because each local error is 
propagated and may either grow or decay as the computation proceeds. 
Decreasing the step size, increases the amount of calculations in the 
program and the accumulated round-off-error can be expected to increase. 
Therefore, as the step size is decreased, to reduce the truncation 
error, the effects of an increasing accumulated round-off-error should 


be considered. 


The subroutine used for the solution of the differential 
equation(s) employed a Hamming's modified predictor-corrector method. 
This method gives an estimate of the local truncation error (d) at each 


step in the program and thus, has the ability to automatically change 
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the step size. The user has to initially specify the step size (h) 

and the maximum allowable truncation error (e). In the subroutine, 

if d>e, h is halved, while if d<e/50, h is doubled. The subroutine 
will not allow the value of h to exceed that initially specified. When 
the program was initiated, it was observed that h is reduced from 

its initial value but, within a few steps the value of h has returned 
to its initial value. Therefore, at this point in the program 

d<e/50 for h/2. For the remainder of the program run the value of h 


is held constant. 


The step size and the maximum truncation error were specified 
as 1/128 and 0.01 respectively. Computations carried out with these 
values gave results that were practically the same as the results 
obtained by using smaller values of h and e. Also, one complete 
program run could be carried out in less than four minutes of computer 


time. 


4.5 Summary 


A set of equations has been developed that will describe 
the behaviour of the electronic impedance of the Impatt diode. These 
equations have been cast into a form that permits them to be accurately 
solved by the simplest computational techniques. Solutions to these 


equations have been obtained using both analog and digital computers. 


The solutions obtained on the analog computer were adequate 
to model steady-state conditions but were not sufficiently stable to 


yield accurate results when simulating transient conditions. Therefore, 
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the final solutions to the equations were obtained on an IBM 360/67 


digital computer. 


The structure of the equations describing the Impatt diode 
facilitates the correspondence between the computed and measured values 
of the electronic impedance. In the equations of the computer model, 
there are several parameters that have to be assigned values. It 
has been found that the more critical parameters have sufficiently 
independent effects upon the calculated electronic impedance, that 


they can be adjusted in a straight forward manner. 
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CHAPTER V 
THERMAL PROPERTIES OF IMPATT DIODES 
5.1 Introduction 


In avalanche diodes, the temperature of the p-n junction 
influences the instantaneous magnitude of the junction breakdown 
voltage and hence the instantaneous magnitude of the junction current. 
Since the magnitude of the diode's electronic impedance is a 
sensitive function of the magnitude of the junction current, the 
characteristics of pulsed, avalanche diodes will change during the 
pulse on-time, as the junction temperature rises. In avalanche diode 
(11) 


oscillators both amplitude and frequency modulation effects occur 


as a result of these temperature fluctuations. 


Previously, the junction temperature characteristics of 
these diodes have been poperted ume) ts within a few hundred nanoseconds 
of the initiation of the bias voltage pulse. In Digital Communications 
pulse durations of tens of nanoseconds may be used. Consequently, 
much more detailed information will be required on the time dependence 
of the junction temperature. The temperature-sensitive property of 
the semiconductor junction that is most directly related to the 
operation of the Impatt diode is the reverse-bias breakdown voltage. 
This temperature-sensitive property can be utilized to obtain a re- 
lationship between the junction temperature and the breakdown voltage. 


Then, if the time dependence of the breakdown voltage can be measured, 
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the time dependence of the junction temperature may be determined. 
The relationship between the junction temperature and breakdown volt- 
age is determined with steady-state measurements, while the time 
dependence of the junction breakdown voltage requires measurements 


taken under transient conditions. 


9.2 Steady-State Thermal Properties 


As the junction temperature increases, the carriers 
interact more strongly with the crystal lattice and require a greater 
field strength to gain sufficient energy for impact ionization to occur. 
Therefore, the junction breakdown voltage increases with increasing 
temperature. The relationship between the diode junction temperature 
and the breakdown voltage was determined by measuring the breakdown 
voltage as the packaged diode was heated in an oven. Measurements 
were recorded for constant temperature levels ranging from 20°C to 
150°C using the measurement system shown in Fig. 5.1. During the measure- 
ments, a constant direct current of 100 microamperes was maintained 
in the diode. This was necessary to ensure that microplasma effects 
were minimized and avalanche breakdown occurred across the entire 
junction area, thus ensuring that the junction was uniformly heated. 
This level of direct current results in only about 10 milliwatts of 
power dissipation within the diode junction and is small enough to have 


only minor effects on the measurements. 


The effects of non-zero power dissipation in the junction 
can be accounted for if the thermal resistance between the diode junction 


and the heat sink is known. As shown in (32), the thermal resistance 
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(r) is defined as 


Tener (1) 
where Tin is the average junction temperature associated with the 
average junction power dissipation Le a Then the corrected junction 
temperature Moe" is given by 

aes = ihe 7Y Vg(75)1 (2) 
where Vg (7) is the measured breakdown voltage corresponding to the 
heat sink temperature 1; and I is the direct current at the time 
of the measurement. Thus, results of measurements of breakdown volt- 
age versus temperature can be corrected for power dissipation within 


the junction itself. 


The junction temperature-breakdown voltage measurements 
were carried out using commercially available Hewlett Packard 5082-0437 
Impatt diodes. Results obtained for two diodes were practically 
identical and the typical measured temperature dependence of the 
breakdown voltage is shown in Fig. 5.2. From these measurements r 
was found to have a numerical value of 22.5°C/watt and the slope of 
the Vg (75) versus ee plot was linear from 20°C to 150°C, with an 
average slope of 0.12 V/°C. The results obtained in these measurements 
are in close agreement with those of previous work 32) The pertinent 
characteristics of the diodes used in these measurements are given in 


Table 5.1, where the values were taken at room temperature and the 


depletion region width and doping levels were derived from 1 MHz 
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capacitance-voltage measurements. It was also noted from these 
measurements, that the diodes were closely approaching a punch through 
condition when the reverse bias voltage reached the breakdown value. 
Since both diodes displayed practically identical temperature 
dependence of the breakdown voltage, it would appear that differences 
in the depletion region width and doping concentrations of the range 
shown in Table 5.1 do not alter the steady-state thermal character- 


istics of these diodes. 


Waele Sig 


CHARACTERISTICS OF IMPATT DIODE 
TYPE H.P., 5082-0437 


DIODE} BREAKDOWN DEPLETION DEPLETION AVERAGE 
VOLTAGE LAYER LAYER DOPING 
(100ua) CAPACITANCE WIDTH CONCENTRATION 


9.13x10!° 


6.85x10!° 


5.3 Transient Thermal Properties 


Since little is known of the transient response of the 
junction temperature of Impatt diodes, it was deemed advisable to 
study its behaviour experimentally prior to developing the theoretical 
model of heat flow in the junction area of the diode. The philosophy 


(33) 


of the work reported by Mosekilde et ai in their theoretical 


study on the thermal response of pulsed Gunn diodes, may be applicable 
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to Impatt diodes as well and provide an indication of the results to 


be expected. 
5.3.1 Measurements of Transient Heat Flow 


A double-pulse measurement method similar to that presented 
by nig oe was utilized to measure the time dependence of the 
junction breakdown voltage. Then, with reference to Fig. 5.2, the 
time dependence of the junction temperature was determined. This 
double pulse method utilizes a main pulse of appreciable duration to 
reverse bias the Impatt diode beyond breakdown and a sampling pulse 
of much shorter duration and opposite polarity to return the bias 
voltage to its breakdown level. As the junction temperature changes, 
so will the magnitude of the breakdown voltage and hence, the magnitude 
of the sampling pulse required to return the junction to its break- 
down condition. In order to conduct measurements to within a few 
nanoseconds of the onset of the main bias pulse, the sampling pulse 
must have a fast rise time and be of short duration. This presents 
some problems in designing the circuitry so that it will accept the 
fast rise times and fully react within the duration of the sampling 
pulse. The measurement system is shown in Fig. 5.3, where the coupling 
of the main and sampling pulsers, plus the measurement points A and 
B, were incorporated into a pair of airline tees whose characteristic 
impedance was 50 ohms. All interconnecting leads were kept as short 
as possible. (maximum length about 3 inches) and were constructed from 
heavy copper wire (12 gauge). The d.c. power supply was shunted with 
low parasitic inductance capacitors to enhance the system response 


time. 


batis2979 tad od. cal cae 59m 92 TUG: e  duse 
ans to sonebreqsb amis sav sivesom 03 post tau ae 
oft 8.2 .otd- 09 ‘oonenatst dstwe: nant! .spstiov nwobasend nctiomut 
sent ‘hanimaseb 26 Swwtsisqnet nofsonut ond to gonabnaqeb emtt 
o3 nofIewb sidstosiqas to s2fuq fem 6 2ostlrau port2on e2Tuq sTduob 
seine pntfamse 6 bas nwobdssyd bnoysd sbotb sisqml srs 2e6td sz719v8" 
std odt nwisy of yitrefoq stteoqqo bas noftetub ysdrone doum to son 
-aspneno siuiereginss nottonut siz 2A . evel awobdsevd est of spstiov ‘ 
ebutinpsm aft sone bas apstfov nwobassid alt to sbuttapem std rfitw oz 


-"gaqd 23t od nottonut aid myugay ot bettupey seluq entiqnse 8 
Wo? s ntdtiw of 2dnomesvesem toubnod of vsby0 nl aoe 
gefuq patiqmse ats ,setuq eptd atem ait to tseno onz To 
eingesiq 2tdT .notiewwb tyodz to sd bre omtd sett ses? 6 even Jat 
oni tyso06 {tw tf Jedd oe yittuorta sit potrapresh at ame Tdo7q. 
ontfgmse oft to notte wh oft ntdsiw Jossy yilut bas eomts get szst 
oni fquos sit evodw ,£.2 .ora of nwode ef mateye Inomeywe6om ot: me 
bis A etntog tnamew2ssm aft eufq .21s2fuq pnt iqmse bie nism edd to 
stdztistosted2 920nw e99F satiate to vteg 6 tnt bedsroqraont e190 a 
trode 26 Jqg9% stew 2ebso! pnitoannosyesnt TrA .emlo 02 26w sonsbeqmt 
mov? betowitenos ovowbns (2edont-£ 4u0u6 tena! mumtx een) ‘otd¥ea0q' 2a” 
ittw batnude mete aks 14 yest _. (spun St) svtw ws 


EXTERNAL TRIGGER FOR 


DIODE 
HOLDER 
AND HEATER 


SAMPLING 
PULSER 


502 


TRIGGER 
PICKOFF 


DIGITAL 
VOLT 
METER 


MAIN 
PULSER 


502 


FIG. 5.3 COAXIAL SYSTEM FOR MEASUREMENT 
OF TRANSIENT JUNCTION TEMPERATURE 


87 


CHANNEL 
REAL TIME A 
OSCILLOSCOPE 
500 MHz 
BANDWIDTH CHANNEL 
B 


OSCILLOSCOPE AND MAIN PULSER 


Deel 
POWER 
SUPPLY 


LOW 
PARASITIC 

INDUCTANCE 
CAPACITORS 


wOl 
JIT L2ARAS 


TOUGHT 
TOASAD 


3 


88 


Since the double-pulse method, as used here, allows the 
observation of the main and sampling pulses as they are applied to the 
diode (measurement point A in Fig. 5.3), as well as the current 
through the diode (measurement point B in Fig. 5.3), the pertinent 
pulse shapes are known. The circuit was adjusted so that it would 
fully respond to the sampling pulse in less than 6 nanoseconds. There- 
fore, a sampling pulse with a rise time of 1 nanosecond and a duration 
(exclusive of rise and decay times) of 6 nanoseconds was used. The 
main bias pulse had a rise time of 5 nanoseconds and a duration of 
200 nanoseconds. To ensure that the diode junction temperature had 
closely approached that of the heat sink, between main bias pulses, 


the pulse repetition rate was set to about 500 pulses per second. 


The measurement procedure was to adjust the magnitude of 
the main bias voltage pulse until the diode current, during the main 
pulse, was at the desired level. Then, at selected time intervals 
following the onset of the main bias pulse, the magnitude of the 
sampling pulse was carefully adjusted until the net bias voltage across 
the diode, within the duration of the sampling pulse, had just been 
reduced to the breakdown value. This magnitude of the sampling pulse 
and the time interval were recorded and plotted as the change in the 
junction breakdown voltage with time, relative to the onset of the 


main bias pulse. 
The reduction of the net bias voltage across the diode, to 


the breakdown value, was apparent from observing the instantaneous 


diode current. When bias voltage across the diode reaches the breakdown 
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value, the complex current fluctuations associated with microplasma 


preatdowitcd. 


cause the diode current to vary slightly for each pulse 
in the repetitive bias voltage pulse train. The observed effect on 
the oscilloscope trace is a "flutter " instability in the current 
waveform at the position associated with the trailing edge of the 
sampling bias pulse. During the experiments, the magnitude of the 
sampling pulse was increased until this instability in the current 
waveform became apparent. By selecting the magnitude of the sampling 
pulse where the flutter just becomes definite as a reference, it 


was possible to carry out consistent measurements at all selected 


sampling positions within the duration of the main bias pulse. 


Typical waveforms observed during the measurements are 
Shown in Fig. 5.4. In these photographs, the upper trace represents 
the diode current, while the lower trace shows the main and sampling 
bias voltage pulses. Both the diode current and the main bias pulse 
increase in magnitude toward the bottom of the photograph. The 
magnitude of the sampling bias voltage pulse increases toward the 
top of the photograph. Fig. 5.4(a) shows the overall relationship 
between the main bias voltage pulse, sampling bias voltage pulse and 
the diode current response to both these bias pulses. Fig. 5.4(b) 
gives an expanded view of the effect upon the diode current when the 
sampling bias voltage pulse is applied. In Fig. 5.4(c) the magnitude 
of the sampling bias voltage pulse has been increased until the character- 


istic "flutter" in the diode current is apparent. 


jnsiws ont nt cavebaa it : watt one ee : 7 
gid to opba pnttters sft dtiw besefoores nottteog ott 36 motavew © 

ent: Yo abutinpem any ,2inamiisqxs ort pntrud .s2fuq estd ont fqmez ; 
tiisawwo odd nt ysrfidssent ett friqu beesstoaT Baw szfuq brit Feqse 
pat fqme2 atid to sbusrapsnt aif patsosise ya .Insvaqas oms2od egviles 
gt ,Sonoyvetey 6 26 otfnrteb zemiaoed t2ut wttult odd svedw s2fuq” 
bojneise ffs is 2dnomew2eom tnatefeno> Io YG ot efdteeoq 26w 


sefua 2std nfem ons Yo notiswb ahs ntditw anotti2oq on t'Tamnee 
‘4 oe 
sis 2snamsyu2sem edt pniyub bayvisedo emvotevew fsatqyT 


sinpesigey S76713 TEqqu sts .2dgsipotodg sestid nT .8.e .pra at nwone 
prffams2 bre nism ald ewoke aosid Yswol sant sftiw ,inevwo abotb ont 
s2fuq 2atd nism sid brs tnaviws sbokb ott diod .zeeTuq opsdtov zeta 
sit rigsyposoilq) rtd Yo moddod’ snd btewot shud taped nf sesetont 

sft biswot 29285 70nt s2fug spas fov estd ont lamse sfz to abut tngem - 
athenotssloy Tferavo sit awode (5)s.2 .pt7 .riqerpotorg ong 0 got 
bas gzTuq spetfov estd entiqmse .saluq spsdfov 26td atem ony ngewied 
(d)hi@ ptt .zeeluq astd e2edt dad of senoqasy InervWw Ssbotb ott 
of? warkw tasyw> sbotb sit noqu ¢397%9 sit to woty pabnsqks A ie 280he! = 
abudtopan odd (2)$.2 .eFi nt .betiqqs 2h safuq spedfov whiseps 


90 


Fig. 5.4 PHOTOGRAPHS OF TYPICAL OSCILLOSCOPE TRACES OBSERVED 
DURING THE TRANSIENT JUNCTION TEMPERATURE MEASUREMENTS. 


In all photographs, the upper trace represents the d.c. diode 
current resulting from the application of the main and sampling bias 
voltage pulses. These bias voltage pulses are shown in the lower trace. 
The upper trace has a vertical scale of 20 MV/cm and the horizontal 


scale is indicated in the lower right hand corner of each photograph. 


Fig. 5.4(a) Overall relationship between the main and sampling 


bias voltage pulses and the resulting diode current. 


Fig. 5.4(b) Expanded view of the sampling bias pulse and its 


effect on the diode current. 


Fig. 5.4(c) Characteristic flutter in the diode current 
waveform associated with the reduction of the bias voltage to the breakdown 


level. 
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From Fig. 5.4(a), it is observed, that, as time measured 
from the onset of the main bias pulse increases, that the diode current 
gradually decreases and the main bias pulse voltage gradually increases. 
This decrease in the diode current is a result of the increasing 
junction temperature, while the increase in the main bias pulse is a 
Characteristic of the measuring setup. With reference to Fig. 5.3, 
monitor positions A and B were connected to the 50 ohm impedance 
oscilloscope channels A and B, respectively. Thus, the voltage applied 
to the diode, monitored at position A, is the fixed voltage pulse 
supplied by the main pulser minus the voltage drop across 3-50 ohm 
impedances in parallel. As the diode current decreases, so does the 
voltage drop across the 16.67 ohm impedance and the resultant voltage 
at A shows a corresponding increase. This effect may lead to some 
inaccuracy in the measured results; however, this effect should be 


very small. 


The waveforms shown in Figs. 5.4 display some character- 
istics which are not an integral part of the measurements. The small 
overshoot at the leading edge of the main bias pulse is a character- 
istic of the circuitry and the main bias pulser, consequently part of 
the initial rapid decrease in the diode current shown in Fig. 5.4(a) 
is due to the decay of the overshoot in the buildup of the main bias 
DUISG, AISOsuinarigense4(C)), stnesinicialerecurn spike of Current 
coinciding with the trailing edge of the sampling pulse is due to 
capacitive charging in the circuitry. The actual buildup of diode 
current after the decay of the sampling pulse is represented by the 


unstable portion of the waveform. 
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Measurements of the increase in junction breakdown voltages 
with time, were carried out at various heat sink temperatures in the 
range 25°G to: 1507€2 Slo avoid the complex current fluctuations 
associated with microplasma breakdown at the onset of the main bias 
pulse, a diode current of 1 ma was maintained before the application 
of the main bias pulse. Complete tests were conducted at each base 
temperature using two separate Impatt diodes and two magnitudes of 
the main bias pulse. The main bias pulse was adjusted so that average 
diode currents of 15 and 30 ma were maintained for the duration of 
the main bias pulse. Each test, within itself, was very consistent 
with little scatter among the plotted measurement data. The experi- 
mental results showed little change as the heat sink temperature was 
increased from 25°C to 150°C. While the diode current level did decrease 
as the heat sink temperature was raised (for a fixed magnitude of main 
bias pulse), the variation of the diode junction breakdown voltage 
stayed relatively the same. For the same reference conditions, the 
diodes used here produced generally similar results, with only small 


differences in the base levels and slopes of the plotted results. 


Typical behaviour of the junction temperature with time is 
shown in Fig. 5.5. Also shown in this Figure are the theoretical results 
which are included here for comparison purposes in later discussions. 

The two sets of curves correspond to the two current levels used for 
the tests, namely 15 and 30 milliamperes. In Fig. 5.5,the zero time 
corresponds to the onset of the main bias pulse, which is taken as 

the point on the initial voltage buildup at which the magnitude of the 


pulse equals 90% of its maximum value,including the overshoot, Previously 
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reported (22) 


measurements of the temperature dependence of the diode 
junction can be interpreted to within 100-200 nanoseconds of the onset 
of the main bias pulse. These earlier measurements are consistent 
with the measurements reported here, in that the final slopes of the 


plotted measurements given in Fig. 5.5 closely agree with the initial 


slopes of the previous measurements. 
5.3.2 Thermal Model For Transient Heat Flow 


The heat generated in a p'n Impatt diode by the direct 
current power dissipation is largely confined to the high resistivity 
region, with the major part of the heat being generated near the 
junction where the electric field has its largest magnitude. The 
argument used to justify one-dimensional analysis for the elect- 
rical properties of the Impatt diode can also be used to justify one- 
dimensional analysis for the thermal properties of this diode.Confining 
the discussion to one-dimensional analysis, the heat flow in the diode 
can be represented by a semi-infinite solid of length &% with heat 
generated uniformly and at a constant rate within a region at one end 
of the solid. This model is shown in Fig. 5.6, where the heat generat- 
jon is confined to the region o <x¢ y.With reference to Fig. 5.7, which 
Shows the physical structure of the radeon by taking the x=o plane 
as the center of the heat generation region, the boundary conditions 
at this reference position are those of a thermally insulated plane. 
Far from the heat generation region the temperature will be held at 


some constant value, such as that of the heat sink. Under these 


conditions the heat conduction equation is 
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AN H.P. 5082-0437 IMPATT DIODE 
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where the variation in time (t) of the temperature (T) at various 
positions (x) in the solid, depends upon the thermal diffusivity (D), 
the heat capacity (Cy) and the rate of heat generation (H) per unit 


volume, per unit time. The applicable boundary conditions are 


also, the parameters D and Cy have been assumed to be independent of 
temperature and constant throughout the solid. Since only the change 
in temperature with position and time is of interest, the initial 
temperature (T,) may be taken as zero. With these boundary conditions, 


the solution to the heat conduction equation'?>’ given in Eq. (1) is 


T = a 1 - 2 i* erfc [ fy-x) ]- 2 i* erfc p (ytx) jhe soazugeneyt(3) 
V Vv4Dt V4Dt 
and 
T= at, i* erfc [ (cy) 5 - i* erfc px) ] x>y (4) 
V V4Dt V4Dt 
This solution is essentially the same as that given by Mosekilde®9 o¢ al. 
and when taken at X=o can be written as 
T = ae 1-4 i? erfc [ eh ey }} (5) 
V 2/ 2 


This last equation is not bounded for t.+ » and will not accurately 
describe the temperature for steady-state conditions; however, the 
transient heat flow under pulsed bias conditions should be adequately 


represented. 
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These equations will describe the heat flow in either direction 
in the Impatt diode (Fig. 5.7) if the diffusivity constant is similar 
for both directions and if there is a very low thermal resistance 
between the silicon and the copper. The diffusivity constant (35) for 
Silicon and copper are given as 0.8 and 1.14 cm/sec, respectively, 
for room temperature conditions. From the steady-state temperature 
measurements, the value of thermal resistance between the diode junction 
and the heat sink was found to be low (22.2 °C/watt). This thermal 
resistance includes the effects of the bulk semiconductor, and the 
heat-sink-diode-package union plus the interface between the semi- 
conductor chip and the copper stud. Therefore, the thermal resistance 
of this silicon-copper interface will be significantly less than 
22.5°C/watt. Consequently, Eqs. (3), (4) and (5) should adequately 
describe the heat flow in the Impatt diodes used in the measurements 


reported here. 


Since the p-n junction of the diode will be very close to the 
x=o plane of the heat flow model, the time-dependent temperature profile 
of the junction can be represented by Eq. (5), which is plotted in 
Fig. 5.5. In this Figure, the parameter A, which represents the quantity 
H/Cy » was chosen for a close correspondence of the theoretical to the 
experimental results. The experimentally observed initial rise of 
junction temperature is greater than the theory would indicate. This 
could be accounted for by some small delay in the initial transfer of 
heat away from the junction area or some inherent inadequancy in the 


one-dimensional approximation of the heat flow model. Otherwise, 
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there is generally good agreement between the experimental results 
and the theoretical model for the transient heat flow in the junction 


area of the Impatt diode used here. 


When considering the heat flow through regions of the semi- 
conductor diode which lie outside the immediate vicinity of the 
junction, Eqs. (3) and (4) have to be considered. These equations 


have been plotted in Fig. 5.98/33) 


as a function of the parameter 

Dt/y. Considering the physical structure of the Impatt diode given 
in Fig. 5.7, an appropiate value for y would be 2 microns, then using 
D=0.8 cm/sec for silicon and t=200 nanoseconds, the parameter Dt/y* 
is evaluated as 4.0. With reference fo Figs.5.7 and 5.8, it is 
observed, for durations of the main bias pulse of 200 nanoseconds or 
less, that heat flow from the junction area has penetrated only about 
half the n* region. Therefore, the relatively poor heat sinking at 
the wire bond should have little influence on the junction temperature 
for the low repetition rate of pulsed bias operation used here. The 
heat flow through the p region has only a few microns distance to 
travel before reaching the heat sink. Hence, the heat flow through 
the p region will be the governing factor in determining the junction 
temperature. 


From the results given here and those reported earlier! 32>33) 


the simplified thermal model, represented by Eqs. (3), (4) and (5) would 
appear to adequately describe the heat flow in the region of the semi- 
conductor junction under pulsed bias conditions. It is likely that the 


thermal model only describes the heat flow during the first 30-40 nano- 
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seconds and then the heat flow from the relatively small silicon-copper 
interface, to the larger copper stud, becomes increasingly important. 
However, this iS a more complex problem where more than one-dimensional 
analysis may have to be used and is beyond the scope of the work 


being carried out here. 


0.4 Computer Modelling of Thermal Properties 


From the steady-state thermal properties of the Impatt diode 
it is known that the junction temperature is directly related to 
the junction breakdown voltage. Thus, for a given value of reverse 
bias voltage across the diode junction, an increase in the junction 
temperature of a certain amount is taken to correspond to a known 
decrease in the effective bias voltage applied to the junction, 
therefore, the effect of junction temperature changes can be incorporated 
in the equations of the computer model as a change in the d.c. bias 
voltage applied to the diode. With reference to the normalized 
equations given in the previous chapter, Eqs. (7) and (14) Section 
4.2, the term Vp in the electric field equation represents the 
bias voltage applied to the diode and it would be this term that 
would be varied with time to account for the changes in the junction 


temperature. 


To provide an accurate comparison between the computed and 
experimental results, both the decreasing current through the diode 
and the increasing main voltage pulse across the diode (reference Fig. 
5.4(a)) were taken into account. In the computer model the bias voltage 


applied to the diode junction (Vp) consisted of three components, the 
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fixed bias voltage set at the main pulser (V_)> the voltage drop 
across the 16.67 ohm resistance, at measurement point A, due to the 
current flow in the circuit (Vp) and the change in the breakdown 


voltage due to the increasing junction temperature (V The net 


1): 
bias voltage applied to the junction is denoted Vp and can be written 
as 


Ve =v =p = Ve (1) 


where the values of Vy are taken from the experimental results and 


both Ve and Vy will change with time. 


In arriving at a relation that dont describe the time (t) 
dependence of Ves Eq. (5) sub-section 5.3.2, was replaced with a 
closer approximation to the experimental results by using a vt 
dependence from 0 - 40 nanoseconds and a linear dependence on t from 
40 - 200 nanoseconds. Following the same procedures as outlined in 
Chapter IV, where the scaled time (t) was used, the following relations 


for Vy were obtained: 


- for o < t < 40 nanoseconds o< ts 400 
Vz = 0.0057 Vt for the 15 ma current level 
V, = 0.0084 VT for the 30 ma current level 

- for t > 40 nanoseconds t > 400 
Vy ='0. 00011 for the 15 ma current level 
Vy = O200017'¢ for the 30 ma current level 


When these values were applied as Le the computed results were found 
to be practically identical to the experimental results in the time 


variation of both the diode current and the main bias pulse (reference 
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Fig. 5.4(a)). The computations were carried out using room temperature 


data for the ionization rates. 


9.5 Summary 


In this chapter, both the steady-state and the transient 
thermal properties for a type of Impatt diode have been investigated. 
A double-pulse method has been used which permitted measurements of the 
junction temperature to be carried out to within 5 nanoseconds of the 
onset of the bias voltage pulse. These studies have related the 
experimental measurements to a theoretical model for heat flow in a 
semiconductor and also to the computer modelling of the Impatt diode. 
Good agreement has been obtained between these three aspects of this 


study. 


A study of the time dependence of the junction temperature 
was carried out by measuring the time dependence of the junction 
breakdown voltage and the steady-state relationship between the junction 
temperature and breakdown voltage. These two sets of measurements 
were conducted over a range of temperatures from 20°C to 150°C. It 
was found that the relative change in junction temperature, during the 
application of a bias voltage pulse, was essentially the same for the 
entire range of heat sink temperatures used here. By extending the 
measurements to include the time interval of 100-200 nanoseconds after 
the onset of the bias pulse, it was possible to compare the results 
with previously published workece? The resulting comparison showed 
good agreement in the change of junction temperature with time. A 


temperature of 150°C was the upper limit of the measurement system as 
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it was used here. The main limiting factor is the method used to hold 
and heat the packaged diode. Electrical interconnections have to be 
kept short to preserve the system response to the fast rise time pulses, 
while at the same time, the high heat sink temperatures should not 

be applied to the bias voltage pulsers. Extensive design work may 

have to be carried out if heat sink temperatures greater than about 


150°C are to be used with this method. 


A theoretical model of the time dependence of the junction 
temperature was arrived at by considering a one-dimensional analysis 
of a semi-infinite semiconductor, with a well-defined location of heat 
generation at one end. This heat flow model gave good agreement with 
the experimental results but the model will be limited to representing 
Only pulsed operation of the Impatt diode. It has been found, for the 
low pulse repetition frequencies used here, that the poor thermal 
path presented by the bonding wire in the diode package will have minimal 
effect on the junction temperature. However, there will be an upper 
limit on the pulse width and repetition rate, above which the thermal 
path presented by the bonding wire will affect the time dependence of 
the junction temperature. Then, a much more detailed analysis of the 


heat flow in the diode and its package would have to be carried out. 


Incorporating the time dependence of the junction temperature 
into the computer model of the Impatt diode was readily carried out. 
A simple variation of the simulated d.c. bias voltage sufficed to 
reproduce the experimentally observed variations in diode current and 


bias voltage. Since the junction temperature mainly alters the diode 
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current, as long as the computed and experimental results are compared 
On the basis of similar diode currents the effects of junction temp- 


erature changes on this comparison, should be minimized. 


In the following work on transient impedance studies of 
Impatt diodes, a knowledge of the behaviour of junction temperature 
as a function of the duration of the bias voltage pulse will be 
important. The buildup of R.F. oscillations in a diode-coupled cavity 
can be of the order of 100 nanoseconds, in which time the junction 
temperature of the diode may change appreciably, leading to changes in 
the electronic impedance of the diode and thus frequency changes in 
the R.F. If the bias pulses are of short duration (less than 40 
nanoseconds) then the effects of the junction temperature change will 
be significantly reduced and may even be able to be ignored in 


certain circumstances. 
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CHAPTER VI 
TRANSIENT MICROWAVE IMPEDANCE MEASUREMENTS 
6.1 Introduction 


Impedance measurements provide an important means of studying 
the microwave properties of Impatt diodes and are essential when designing 
amplifier and oscillator circuits for these diodes. Several authors (3/740225) 
have published methods and/or results of microwave impedance measurements 
On packaged diodes. All of these measurements have been made under 
steady-state conditions. The usual procedure is to measure the admittance/ 
impedance of the packaged diode at some convenient reference plane, and 
then transform the results through an appropiate equivalent circuit to 
arrive at the admittance/impedance of the semiconductor diode chip itself. 
The impedance measurements mentioned above are all performed by measuring 
the Voltage Standing Wave Ratio (VSWR) or the voltage reflection coefficient 


(r) in a transmission line . (The general relations between VSWR, rT and 


impedance are given in Appendix B). 


Even though a lot of work has been carried out in improving the 
microwave measurement techniques for Impatt diodes, they are still limited 
in their absolute accuracy. The magnitude of the packaged diode's negative 
resistance is of the order of 10 ohms, or even much less, for practical 
diodes. Under conditions of low d.c. currents, high R.F. voltages or 
maximum power output, this magnitude can be as low as 1 or 2 ohms. 
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VSWR or reflection coefficient measurement techniques are sensitive 
enough to produce highly accurate results at these low magnitudes. They 
advocate the use of a microwave impedance bridge technique that requires 


a very stable frequency source and some critically machined components. 


The accuracy of the measurements is further restricted because 
the Impatt diode is a non-linear device and some harmonics of the funda- 
mental frequency will be produced during the active measurements. Most 
measurement equipment does not possess sufficient bandwidth to handle the 
second harmonic component. Also, even small discontinuities in the 
measurement equipment and transmission lines used can cause the second 
harmonic component of the R.F. frequency to become significant at R.F. 
voltage magnitudes of only a few volts (25) Therefore, all large-signal 


microwave impedance measurements were performed with R.F. voltage levels 


in the range of one volt (r.m.s.). 


To date, a detailed comparison on the accuracy of the various 
applicable impedance measurement methods has not been carried out, to 
this writer's knowledge. However, the standard measurement methods have 
been used, with satisfactory accuracy, to perform similar impedance measure- 
ments to those required for the Impatt diode. Therefore, it is reasonable 
to expect that presently available measurement methods are sufficiently 
accurate to uncover the important behavioral characteristics of the diode's 


microwave impedance. 


The previously-mentioned measurement methods are suitable only 


for steady-state operating conditions. Under transient conditions, the 
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diode impedance is a function of time, and as a result, the measurement 
techniques used must have the capability of displaying this time-dependent 


impedance over the time and impedance range of interest. 


6.2 Transient Impedance by Direct Comparison of Voltage/Current Waveforms 


This method (42»43) 


requires recording of the R.F. voltage wave- 
form across the diode as well as the current waveform through the diode. 
These recorded waveforms are then Fourier analyzed and the fundamental 
frequency components compared to yield the desired admittance or 
impedance. If these waveforms can be monitored during the rise and fall 
times of the bias voltage pulse, the transient behavior of the diode's 
impedance can be visually observed. This method would be very useful for 


monitoring the build up of oscillations in a cavity and the effect of 


this build up of R.F. oscillations on the coupled Impatt diode. 


In theory, this general type of impedance measurement method 
is very versatile. However, its usefulness is restricted by the practical 
difficulties that arise when the method is actually implemented. The 
"cavity oscillator method" and the “airline tee method" to be described 
below, are discussed to emphasize the practical problems of conducting 
transient impedance measurements and also, to illustrate the limitations 
of actual measurements that are made. Both of these methods have been 


tested and their limitations, along with suggested solutions, are presented. 


6.2.1 Cavity Oscillator Method 


This particular measurement method is very versatile, as it is 
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possible to carry out impedance measurements under a variety of operating 
conditions. In this method the diode is placed in a cavity tuned to a 
specific frequency. Then, the R.F. voltage in the cavity and the R.F. 

current through the diode are monitored on a sampling oscilloscope 
alternately, an oscilloscope-driven chart recorder can be used to give an 
enlarged presentation of the waveforms. The waveforms are then directly 
compared in magnitude and phase, or Fourier analyzed, to give the admittance/ 
impedance measured. The use of a simple, fixed frequency, cavity enables 

the measurement set up to be utilized in either the R.F. driven state 


or the self-oscillatory state. 


In the REF. driven state, the cavity is sufficiently loaded 
so that the coupled diode will not self-oscillate and R.F. voltage is 
supplied by an external source, via the part which normally is the cavity 
output. When the self oscillatory state is employed, the cavity is un- 
loaded to the extent that self-oscillations will occur. Then the loaded 
Q of the cavity can be varied and the build up of oscillations observed 


for the various conditions of the external circuitry. 


Utilizing a combination of these two operating states, extensive 
experimental measurements can be made. The measurement system is shown 
in Fig. 6.1 and an expanded view of the cavity is given in Fig. 6.2. To 
ensure that only the TEM mode exists in the cavity, even at the second 
harmonic of the R.F. frequency, 7 mm airline was used for the cavity 
construction. These diode cavities were of the single-transformer, 


fixed-frequency type, so that they could be represented by a simple 
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equivalent circuit. A change of operating frequency required a change of 
center conductor (due to spring loading the diode contact); however this 


was a straight forward procedure. 


The R.F. current waveform through the diode was monitored by 
the voltage developed across a one ohm disc resistor placed in series with 
the packaged diode and the copper heat sink. Monitoring of the R.F. 
voltage was carried out by a capacitive probe located in the vicinity of 
the diode. The system could be calibrated by measuring the packaged diode 
impedance at selected bias voltages (below and just at breakdown) in 
the cavity and then comparing these results with carefully made VSWR 


measurements of the packaged diode by itself. 


Although this measurement method has a wide scope of possible 
uses, it suffers from practical limitations. The more restricting of 
these limitations are outlined below: 

— For accurate measurements, it is necessary to place the 
current viewing resistor as close as possible to the 
packaged diode. However, locating the resistor between 
the packaged diode and the heat sink places an upper limit 
on the steady-state current level through the diode. This 
is due to the relatively high thermal resistance through 
the ceramic base of the disc resistor. It is not practical 
to reverse the diode package, as then the center con- 
ductor of the cavity would need to act as the heat sink 


and its thermal resistance is not satisfactory either. The 
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upper limit on the steady-state current flow is there- 
fore only a few milliamperes which, for most diodes, is 
not sufficient for self oscillation to occur. By using 
special ceramic bases for the thin film resistors, such 

as Beryllia , the thermal conductance can be improved, but 
is still lower than would be the case when the disc is 

not present. Thus, the d.c. bias current is limited to 
values sufficient to produce low-level C.W. oscillations 
and, in the pulsed mode, the diode is limited to low duty 


cycles. 


— To provide sufficient sensitivity, the depth of penetration 
and the diameter of the voltage probe have certain minimum 
values. While to ensure that the entry port of the probe 
does not unduly disturb the electric field in the cavity, 
this entry port has a maximum diameter. The result is 
that the coaxial impedance presented to the cavity by the 
entry port of the voltage probe and the physical presence 
of the probe itself, constitute sufficient loading of 
the cavity that self-oscillations are not possible. There- 
fore, the measurement system is limited to the R.F. driven 


mode of operation. 


* 
Al Si Mag Beryllia Ceramics (trade name) America Lava Corporation 


Chattanooga Tennessee. 


eT, es s20n 10% ,dotr 
pnteu ye .wo90 ot 5 eel Rio a2 ot gnototttue 4 
Hove jevoreteay mite nba sitt vot esse 3tms199 ke 


tid ,bavoigné ed ‘nea Sanetsbbnos ‘Teined) odd . aE Tyree Sa 


2t seth srt Mod sees oiit od bfuew nod? rowol MTtte 2f 
og bewximtT et gneviia e6td .o.b ofS ,2udT . Jase ton 
anoisef{tozo .W.d fevef-wol soubotg ot snstotttwe 2autsv 
yitb Wor of badimrl zt abotb Sit , aban beztuq oft nf .bns 
| .esfoyo 


notteytansg to diqsb add .yitvittense gostatttwe ebivorg oT 
muntrtm ntsivaa svet sdovg spetfov ont to tadsmarh six brs 
edow si3 to voy viine off teNs atuans of ofinW .zsuTsv 
«Mi fvso eft nt blot? ofvtoats ond dwterb yfubmu ton esob 
ét dfu2sv SAT .19tam6ib mumixem § 26d tyvoq yvIne etd 

oft yd ystveo sdt of batnsestq sansbsqmt Isfxeed ott tent 
gonsesig Isofeyiq sid bas sdovg spsdfov sit to tyvog yrsne 

. Yo entbsol snatoittue stustienos ,tlsest adovg sft Yo 
-9"adT .ofdt2zog fon o16 2notteltfozo-7fse ted? yttveo sn? 
navi4b 1A ons oF batimeéh 2f mateye Inomeuwz5sm 9d .o70t 
-notsie1sqgo to sbom 


0teF aes x 7 


petveg ow on 
‘eo wariae ’ 
ye ims 


. -- AD 


awl pat “— 


114 


— In the R.F. driven mode of operation, with low duty cycle 
bias voltage pulses applied, certain problems arise. In 
Order to ensure that the diode will nearly return to the 
nominal heat sink temperature after application of each 
bias voltage pulse, the pulse repetition frequency has an 
upper limit of about 100 KHz. The frequency of the driving 
R.F. is about 5 GHz and this must be "counted down" to 
100 KHz to synchronize the application of the bias voltage 
pulse and the driving R.F. This synchronization is 
necessary so that the display of the R.F. voltage and current 
waveforms can be correlated to the waveshape of the bias 
voltage pulse. It is difficult to accomplish the frequency 
division, or count down, of the driving R.F. voltage, 
from 5 GHz to 100 KHz and still maintain a sufficiently 


high degree of stability to ensure proper synchronization. 


This measurement system is readily set up, with a minimum of 
equipment, and is versatile in nature. However, the inherent practical 
limitations encountered in the actual implementation of the system, restrict 
its usefulness. The difficulty of counting down the frequency of the R.F. 


4 can probably be resolved with commercially 


voltage by a factor of 5 X 10 
available equipment, but the loading effects of the capacitive voltage 


probe will require extensive design work being carried out. 


6.2.2 Airline Tee Method 


This method is similar in principle to the previously discussed 


cavity oscillator system. The R.F. voltage across the diode is compared 
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with the R.F. current through the packaged diode to arrive at the 
microwave impedance of the packaged diode. Although the airline tee is 
limited to the R.F. driven mode of operation, it has the advantage that 
the R.F. voltage across the diode is single frequency and free of second 
harmonics. This is accomplished by terminating the incident R.F. voltage 
in precision connectors and terminations that are readily available with 


the required frequency bandwidth. 


The physical equipment layout is the same as that given in 
Fig. 6.1, only the diode holder has changed. An expanded view of the 
diode holder is given in Fig. 6.3. This diode holding unit is similar 


to that used by Sonne tonic) 


in the measurement of the switching times of 
p-i-n diodes. In the diode holder of Fig. 6.3, the R.F. voltage at the 
position of the packaged diode is monitored by a capacitive probe and 

the R.F. current through the diode is monitored by a current viewing 
resistor. The current viewing resistor is commercially available and 

is specified to have an impedance of one ohm real part, with negligable 
imaginary part, to frequencies beyond 12 GHz. The use of 14 mm airline 
was necessitated by the physical dimensions of the packaged diode and the 
need for a spring-loaded contact in the center conductor of the airline. 
The physical size brane voltage probe and its entry port into the airline 


were the same as that le in the cavity oscillator method. Therefore, 


minimum distortion of the electric field was ensured. 


| 


Similar to the cavity oscillator system, the use of the current 


viewing resistor presents heat sinking problems and the airline tee 
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method is limited to low duty cycle, pulsed operation. However, this 

is not a severe restriction, as transient behaviour can still be studied. 
The trigger synchronizing problems inherent in this type of operation 
can be circumvented by using the modified equipment set up shown in 

Fig. 6.4. Here, the monitored R.F. voltage and current for the packaged 
diode are down converted to about 500 MHz and displayed on a real time 
oscilloscope with sufficient frequency bandwidth. The various delays 

in the measuring circuits can be adjusted so that a simultaneous 
oscilloscope display of the bias voltage pulse and the down converted 
voltage or current waveform is possible. One limitation inherent in the 
use of the down conversion of frequency, is that the second harmonic 


component in the current waveform may be lost. 


During the calibration of the measuring system, it soon became 
apparent that there were serious practical limitations that restricted 
the accuracy of any actual impedance measurements. The more serious of 


these limitations were the following: 


— Due to the necessity of supplying the Impatt diode with 
a high d.c. voltage for reverse biasing, the termination 
end of the diode holder has to contain a d.c. block. 
This means that a small, but finite, VSWR (v1.1) is 
presented to the incident R.F. driving voltage and any 
second harmonic component generated by the diode. A 
similar situation exists at the driving termination due 
to the bias tee. Thus, there is some distortion of the 
electric field across the packaged diode caused by re- 


flections and some component of the second harmonic may be 
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present. This effect could be minimized by careful 


design of the terminating end of the diode holder unit. 


— The combination of the diode holder in the airline center 
conductor and the physical length of the packaged diode 
has the electrical effect of a shunting inductive post 
across the coaxial airline. This post effect is so 
large at the fundamental frequency of interest (~ 5 GHz) 
that it masks the behaviourof the diode during the 
impedance measurements. This inductive post effect can 
be tuned out by placing an adjustable R.F. short at the 
termination end of the airline tee. After the inductive 
post effect was tuned out, the variations in the impedance 
of the packaged diode, with changing bias voltages, were 
readily observed in this measurement system. However, 
it was difficult to arrive at calibration factors that 
would provide a close comparison between the measured 
impedances and those known values for the same bias volt- 
ages. In general, attempts at calibration proved to be 


difficult and unreliable. 


Similar to the cavity oscillator measuring system, the airline 
tee system is restricted in its usefulness by the practical difficulties 
of actually implementing the measuring system. A possible solution to 
the inductive post effect would be to design the center conductor of the 


coaxial airline in such a way that the physical separation between the 
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inner and outer coaxial conductors, would be equal to the height of the 
diode at the measurement plane of reference. This could readily be 
accomplished by using tapers or three-section, quarter-wave-transformers 
to give the necessary impedance matching over the frequency range of 


interest. 


6.3 Transient Impedances by Carrier Reinsertion 


The major advantage of this measurement method is the inherent 
Simplicity of its operation. Although the method is restricted to the 
R.F. driven mode of operation, it has few other limitations. The 
equipment layout is given in Fig. 6.5 and an expanded view of the diode 


holder is shown in Fig. 6.6. 


The general principle of operation is that the packaged diode 
is mounted in a coaxial holder which is attached to a coaxial four-port 
hybrid. This hybrid separates the incident and reflected R.F. voltage 
waves. The reflected voltage wave and a portion of the incident voltage 
wave are separately routed, through switches, to a waveguide four-port 
hybrid (magic tee). In the magic tee the R.F. voltage waveforms are 
applied to a crystal detector. The switching circuitry allows the re- 
flected and incident signals to be applied to the detector either sing- 
ularly or as a sum. The R.F. crystal detector gives a d.c. output volt- 
age whose level is proportional to the R.F. signal sampled by the 
detector. These detected levels are displayed on a real time oscilloscope 
and then recorded with a polaroid camera. The resulting oscilloscope 
traces on the film are then analyzed to give the complex reflection 


coefficient. 
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Thus, the major limitations of the previously discussed 
measurement methods are overcome. The driving R.F. voltage is supplied 
from a free running source so that the problems of synchronizing the 
driving R.F. source and the bias voltage pulser do not exist. The 
current viewing resistor is only used as a means of determining the 
current pulse through the diode to serve as a reference level for the 
comparison of the computed and experimental results. During the actual 
impedance measurements the diode package is in direct contact with the 
copper heat sink. Measurement sensitivity is good and the calibration 
procedure is straight forward and accurate. There are some restrictions 
on the measurement system,as it was used here, but the results obtained 
should have similar accuracy to those taken by VSWR/reflection coefficient 


measurements under steady-state operating conditions. 


The practical limitations of equipment availability and the 
desirability of the coaxial mount for the packaged diode, necessitated 


the use of both coaxial and waveguide components in the measuring system. 
The equipment layout can be divided into two major divisions, one 


consisting of coaxial components and one consisting of waveguide components. 
The microwave aspects of the measuring system can similarly be divided 
along divisions of coaxial and waveguide components and will be discussed 


accordingly. 


6.3.1 Biasing and Separation of Incident and Reflected Signals 


With reference to Fig. 6.5, the biasing and separation functions 


are carried out in the coaxial part of the measurement system. The 
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incident R.F. voltage, plus the d.c. and pulse bias voltages, are supplied 
to the packaged diode by a bias tee. The coaxial directional coupler has 
very little effect upon the system operation and the d.c. open - R.F. 
termination on port-four of the circulator is accomplished by using a co- 
axial to waveguide adapter with a waveguide termination. There are 

twe critical aspects of the equipment arrangement which will be discussed; 
One is the rise time of the bias voltage pulse at the diode holder and 

the other concerns the frequency bandwidth and R.F. isolation of the 


four-port hybrid. 


Given a voltage pulser with sufficiently fast rise time, the 
rise time of the bias voltage pulse at the diode holder is governed by 
the rise time of the associated circuitry. Since the diode holder is 
mated to a conventional 50 ohm APC-7 connector, the measurement of the 
bias voltage pulse at the diode holder can be obtained by directly connect- 
ing the mating APC-7 connector to the 50 ohm input of the oscilloscope. 
By monitoring the direct output of the pulser, (position A Fig. 6.5) it 
was found that, placing the 50 ohm impedance of the oscilloscope 
channel or a conducting Impatt diode, at the diode holder, very similar 
loading effects were observed on the pulser output. Thus, the bias 
pulse voltage measured by the 50 ohm impedance of the oscilloscope was 
a close approximation to the actual bias pulse across the packaged diode. 
Normally, the rise time of the bias voltage pulse at the diode holder is 
about 10 nanoseconds. This rise time can be decreased to about 5 nano- 


seconds by a judicious choice of capacitors placed across the output 
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terminals of the d.c. power supply. This rise time is maintained as 

long as the rise time of the bias voltage pulser is greater than 5 nanosec- 
onds. The bias voltage rise time at the diode holder can be decreased 

to approximately that of the pulser by redesigning the diode holder so 

that the heat sink is electrically insulated from the system ground. 

Then, the bias voltage pulser is applied between the diode and the system 
ground. This was not done in the measurements carried out here, but 


was tested during the work on the cavity oscillator measurement method. 


The characteristics of the circuitry at the microwave fre- 
quencies used are mainly determined by the frequency bandwidth and R.F. 
isolation of the coaxial four-port circulator. As the diode holder is 
constructed using a precision 50 ohm connector, it has excellent 
characteristics (VSWR < 1.05) to frequencies beyond 18 GHz. However, 
the four-port circulator is limited to one octave of frequency bandwidth 
and outside this frequency range it presents a highly reactive impedance. 
Thus, any second harmonic components generated by the Impatt diode 
would be largely reflected from the circulator and confined within the diode 
mount. In this situation, the R.F. voltage incident upon the Impatt 
diode will have some second harmonic component present. This is a 
similar situation to that encountered in steady-state impedance measure- 
ments and will limit the magnitude of the incident R.F. voltage to the 
order of one volt so that any effect from second harmonic components 


will be minimized. 


Another important aspect of the characteristics of the four- 


port circulator is the isolation between ports 1 and 3(see Fig.6.5). This 
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isolation describes the degree by which the incident R.F. is prevented 
from gaining direct access to the path of the reflected signal. The 
isolation of the four-port used here was greater than 20 dB between ports 
1 and 3, a larger degree of isolation would be desirable however, 20 dB 


was found to be acceptable. 


The measuring characteristics of this circuitry at microwave 
frequencies could be improved by increasing the frequency bandwidth and 
R.F. isolation of the four-port hybrid used to separate the incident and 
reflected R.F. signals. This will not be a simple matter, as coaxial 
components will give two-octaves of bandwidth but the isolation will be 
marginal. Waveguide components will give the high degree of isolation 
required but the frequency bandwidth will be narrow. A possible 
improvement may be to use a broad-band, 3-db, directional coupler with 
very high directivity, although this will result in the loss of one 


half of the reflected signal. 


Deoecehal es DeLeCLION subsystem 


With reference to Fig. 6.5, the R.F. detection takes place in 
the waveguide part of the measuring system. The advantages of using 
waveguide components are that precision components are readily available 
and these components have very low R.F. losses.. The most important part 
of the waveguide portion of the measuring system is the four-port 
hybrid (magic-tee) R.F. detector arrangement shown in Fig. 6.7. The 
magic-tee was not specially compensated, but, was tested and found to have 
equal power division from both the "E" and "H" plane ports, while main- 


taining a high degree of isolation between these ports. Only one of the 
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sum and difference ports was used, as it was found that the outputs of 
these ports are very sensitive to impedance mismatches at their opposing 
member. Therefore, one of these ports was match terminated with a high 


quality waveguide load and the other was used for the R.F. detector mount. 


The R.F. detector used was an integral part of a waveguide to 
coaxial adapter and presented a very good impedance match to the magic 
tee (At the measuring frequency the R.F. detector mount had a VSWR 1.05); 
this detector used a conventional broad-band diode. The sensitivity 
of the R.F. detector was increased by applying a small forward bias to 
the detector diode. This increased sensitivity is essentially due to 
improved impedance matching between the detector output and its load. 

The R.F. detector is designed to be loaded by a high impedance (specified 
as > 75 Kilohms) and if a FET probe is used (input impedance of 1 megohm 
Shunted by a capacitance of 3 pf) high signal levels are recorded from 
the detector output. But the response time of the detection sub-system 
is ~ 14 nanoseconds even with detector biasing. This slow response time 
is due to the high R.C. time constant of the detection system. To 

reduce the time constant, the detector was loaded by the 50 ohm input 
impedance of the oscilloscope channel and then biased for best sensitivity. 
With this arrangement, the response time was better than 5 nanoseconds 
and the sensitivity was sufficient. Detector biasing was accomplished, 
as shown in Appendix C, through a BNC tee to avoid any unnecessary 


introduction of additional capacitance. 


The sensitivity of the R.F. detection part of the impedance 


measuring system could be enhanced by using a more sensitive detector 
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diode. R.F. detection over a broad frequency band is not necessary 

as the frequency bandwidth of the magic tee will limit the R.F. signals 
to their fundamental component only. Resolution of the display on 

the oscilloscope and the subsequent pictures was not a limiting factor 
here. The oscilloscope used had a frequency bandwidth of 500 MHz and 
the same channel was used to calibrate the measuring system, and record 
the reference levels for all the detector outputs. The accuracy of 
transferring the data from the polaroid pictures is quite good if the 


transfer is carried out carefully. 


6.3.3 Calibration of Measuring System 


Before any impedance measurements can be carried out, the 
System as a whole must be calibrated. This is accomplished by taking 
the reference plane for the measurements to be the physical location of 
the large flange on the Impatt diode, as shown in Fig. 6.6. This is 
the most convenient reference plane to take, as it can apply to any 


impedance measurements that are carried out on the diode package. 


Placing the APC-7 electrical short in the APC-7 part of the 
diode holder (i.e. at the reference plane shown in Fig. 6.6) the 
magnitude and phase of the incident and reflected R.F. voltage waves 
V. and Le respectively, are set such that 

Val = V4 

\V. + V. | = 0.0 
This calibration routine accounts for the electrical losses in the paths 
travelled by the incident and reflected signals and also is used to set 


the reference level of the incident R.F. voltage. Once this part of the 
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calibration procedure is carried out for a specific frequency, the 


setting of the attenuators and phase shifters is fixed. 


As shown in Appendix B, to calculate the complex reflection 
coefficient from the information supplied by the oscilloscope traces 
two constants need to be evaluated. The values of these constants 
are arrived at by comparing the impedances measured by the reflection 
coefficient system with those given by carefully made VSWR measurements. 
Since the measurement plane of reference on the diode package is the 
same for both types of measurements, a direct comparison is possible. 
The VSWR measurements were taken on Impatt diodes that were reverse biased 
below breakdown. These diodes were the same ones that were used later 
for the transient impedance measurements. Also, these VSWR measurements 
were repeated to ensure that the high values of VSWR present did not 


lead to inaccurate impedance calculations. 


The general procedure was to select an Impatt diode and carry 
out the VSWR measurements at selected values of reverse bias below 
breakdown. Then,this same Impatt would be transferred to the reflection 
coefficient measurement system and, at the same selected values of 


reverse bias voltage, the impedance would be measured. 


A comparison of the VSWR and reflection coefficient measure- 
ments determined the values of the required constants. Then, the tran- 
sient impedance measurements were carried out before the Impatt diode 


was removed from the reflection coefficient measuring system. 
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6.4 Experimental Measurements 


The hermetically sealed package in which the semiconductor 
diode is mounted, serves to protect the delicate p-n junction and also 
to provide coupling from the diode chip to the external circuit. 
Consequently, the semiconductor diode terminals are not accessible for 
direct impedance measurements. In the process of arriving at a value 
of the semiconductor impedance Vas two other impedances need to be 
considered. These associated impedances are the actual measured impedance 
es and the impedance at the outer reference surface of the diode 
package (Z,). The value of Zn is transformed from its reference plane 
to the reference plane of Ly by an equivalent circuit which is determined 
by the measurement circuitry used. The reference plane for Z, is the 
outer surface of the diode package and the equivalent circuit of this 


package transforms Z_ to Zoe 


p 

The use of the 50 ohm mounting system for the diode package, 
Shown in Fig. 6.6, greatly simplifies the transformation of Zn to Zy 
(37,41,44).Since the impedance transformations involve the use of lumped 
element equivalent circuits to represent distributed effects, the 
simpler the transformations to be represented, the more representative 
the equivalent circuit. The mounting system used here has the dual 


advantage of simplifying the calculations as well as contributing to 


increased accuracy of the impedance transformations. 


The Impatt diodes used in this work were commercially obtained, 
Hewlett-Packard model 5082-0437 diodes, encapsulated in the common pill ~ 


with-two-prong package. “(whose outline is illustrated in Fig. 6.8). 
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All measurements were carried out at frequencies in the range of 5.7 to 
8.3 GHz and hence were close to a common carrier microwave communications 


frequency band. 


6.4.1 Derivation of the Diode Package Equivalent Circuit 


As stated in sub-section 6.3.3, the reference plane selected 
for the impedance measurements was the large flange on the diode package 
shown in Fig. 6.6. This location of the reference plane was selected 
as it was also the reference plane of the electrical short used as a 
reference for both the reflection coefficient and the VSWR measurements. 
The equivalent circuit for representing the transformation of the 
impedance measured at the reference plane (75) to the impedance at the 
outer surface of the diode package (Z,) was arrived at by using a metal 
diode of the same physical dimensions as the real diode package. It 
was determined, by VSWR measurements, that the impedance of this dummy 
metal diode, at the reference plane, could be represented by a lumped 
inductor of 0.6 Nh. This value of inductance was constant over the 


frequency range of 5 to 8 GHz. 


The procedure of using a dummy metal diode to obtain the 
relation between Ly and Zo has been used in previous work 44) where 
reduced height waveguide mounting systems were employed. Although the 
higher order modes necessary to match the TEM mode electric and magnetic 
fields to the metal diode will be somewhat different from those required 
for the real diode package, this difference is considered to have a 
small effect on the impedance transformations. The principle of using 


a single, lumped- inductance, element for the transformation from Z.. 


to Z_ is the same as that used by Getsinger ©) for diodes mounted across 
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a waveguide, in shunt with coaxial conductors and at the end of a 


(45) 


shorted coaxial line. Owen and Cawsey also found that the single 


inductive element was sufficient to transform Zi. to Z 


(37) 


0° 
Getsinger also indicates that, as the outer diameter of 
the diode package approaches that of the coaxial inner conductor, the 
radial modes that require matching to the TEM fields vanish. Then, the 
value of the lumped-inductive element is that of a shorted transmission 
line whose characteristic impedance is formed by the diode package 
outer diameter and whose length is equal to the height of the ceramic 
Sleeve of the diode package. In the package mount used here the 
diameter of the coaxial center conductor is 0.120" and the outer dia- 
meter of the package is 0.080", so very little radial mode matching 
should be required. This is also indicated experimentally, as the 
measured impedance of the metal diode is very nearly that of a short 
circuited section of transmission line whose inner conductor diameter 
and length are those presented by the dummy diode. 


The above results are in good agreement with Owens and 


Gaveeye (>) 


work, who carried out very extensive measurements to character- 
ize a diode package of very similar dimencibies The equivalent circuits 
used by Getsinger, Owen and Cawsey, along with the version derived for 

this work are shown in Fig. 6.9. Owens, introduced a capacitance (C,) 

that was due to the ceramic cylinder in the diode package but found 

that this only had significance at frequencies above 12 GHz. If the 


value of C, is introduced at frequencies aroung 6 GHz, its only effect 


is to slightly reduce the magnitude of Zi before it is transformed by 
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/—— DIODE PACKAGE 
f L, IS INDUCTANCE OF METAL POST 


(a) 


REFERENCE PLANE 
FOR 
IMPEDANCE MEASUREMENTS 


L, = 0.64 nh 
Ly = 0.23 nh 
C, = 0.04 pf 
Co SAGEM Ue Tag 
Ly = 0.6 nh 
Lo = 0.4 nh 
C, = OsL5 spi 


t—we DIODE PACKAGE 
(c) 


Bal Gr EQUIVALENT CIRCUIT FOR PACKAGED IMPATT DIODE 


6.9 

(a) GETSINGER'S MODEL 

(b) OWENS AND CAWSEY'S MODEL 
(c) MODEL USED FOR THIS WORK 
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the transmission line inductance L, (L, this work) to Z,- Therefore, 


] c 
the value of C, was omitted for this work. 


The values of Cy and L, (Figure 6.9 (c)) are essentially 
those listed on the diode specification sheet as the package capacitance 
and inductance, respectively. In the diode package shown in Fig. 6.8, 

Cy represents the capacitance between the top of the post and the top 

cap of the package, while Lis mainly the inductance of the lead wire(s). 
The value of Ly can be determined from the difference in the inductive 
impedance of a metal diode and that of an internally shorted diode 
package, when measured at the same impedance reference plane. The 
internally shorted diode package that was available had the same outer 
physical dimensions as the actual diode package but the internal structure 
was unknown and probably slightly different than that of the actual 

diode used. However, these measurements were carried out to obtain a 
representative value for Lo: A value for Ly of 0.35 - 0.4 nh was found 

in the frequency range of 5 to 8 GHz. The manufacturer's listed values 

of Ly and Cy were 0.5 nh and 0.2 pf, respectively. Thus, the values of 

Lo and Cy for the actual diode's package would be close to the above 
values. Using these magnitudes as a first approximation, the final 

values of Ly and Cy were arrived at by fitting the depletion layer 
capacitances, measured at 1 MHz on a Boonton Bridge, to those calculated 
from a transformation of ZL. to Le using the equivalent circuit given 


in Fig. 6.9(c). This matching was carried out for a range of reverse 


bias voltages below breakdown. 
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By choosing the values of Ly and Co to be 0.4 nh and 0.15 pf 
respectively, the calculated capacitance was within 4% of the measured 
values at 1 MHz for bias voltages of 40 to 90 volts (breakdown voltage 
was 100 volts). This completed the determination of the necessary 


equivalent circuits required to transform L, to Zee 


6.4.2 Experimental Results 


In the experimental setup shown in Fig. 6.5, there are three 
monitor positions, each of which is designed to be connected to a 50 
ohm impedance. The monitor position "A" was provided by a tee and 
this position was continuously loaded by one of the 50 ohm channels 
of the oscilloscope. The continuous loading served two purposes. One 
was to present a reasonable match to the 50 ohm output impedance of 
the bias voltage pulser, as the d.c. power supply was a large impedance 
mismatch; the other purpose served was to provide a reference bias 
voltage level so that the measurements would be reproducible. Monitor 
positions"B" and "C", alternately shared the remaining 50 ohm channel 
of the oscilloscope. Position "B" was a 50 ohm connector in parallel 
with the 1 ohm current-viewing resistor. This was used just prior to 
the calibration of the R.F. detector, to enable the current pulse 
through the diode to be obtained. The current pulse through the diode 
was necessary to serve as a reference for comparison between the theoretical 
calculations and the experimental results. Monitor position "C" jis a 
50 ohm tee that shared the output of the R.F. detector with the bias 
voltage supply for the R.F. detector. This provided a positive d.c. 


offset to the oscilloscope that tended to balance the negative output 
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voltage of the R.F. detector and permitted the use of the most sensitive 


range of the oscilloscope channel. 


It was found that, when the diode package holder was loaded 
by the 50 ohm input impedance of the oscilloscope channel, in lieu of 
an avalanche diode, the waveform observed at the bias pulser 
(monitor A) was very similar to that observed at "A", when the diode 
holder held an avalanching diode. This effect is apparent when traces for 
"A" are compared in Figs.6.10(a) and (b). In this way, the rise time of 
the bias voltage pulse at the location of the Impatt diode was determined. 
Also apparent from this measurement was that the current pulse through 


the diode follows the voltage pulse across the diode very closely(Fig. 6.10a)). 


To avoid instabilities in the monitoring oscilloscope traces 
when repetitive bias pulses are applied to the Impatt diode, a steady 
d.c. current of 1] ma was maintained in the Impatt diode. This ensures 
that the Impatt diode has an identical response for each pulse in the 
repetitive train of bias pulses and the oscilloscope traces are steady. 
This d.c. current level is sufficiently low that the temperature of 
the semiconductor junction is essentially that of the copper heat sink. 
The bias voltage pulses had a duration of about 30 nanoseconds and a 
repetitive frequency of 100 KHz, giving a duty cycle of 0.3%. Thus, 
it was assured that the semiconductor junction temperature closely 


approached that of the heat sink,in between bias voltage pulses. 


The procedure followed for each set of measurements started 
with the determination of the Impatt diode current response for each 
specific value of bias voltage pulse to be used. Next, the current- 


viewing resistor was replaced with a copper disc and the reverse bias 
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Fig. 6.10 PHOTOGRAPHS OF TYPICAL WAVEFORMS OBSERVED ON THE 
OSCILLOSCOPE DURING ELECTRONIC IMPEDANCE 
MEASUREMENTS USING THE CARRIER REINSERTION METHOD. 


In all photographs, the vertical scale for current and 
detected R.F. waveforms is 10 MV/cm. The voltage waveforms have a 
vertical scale of 500 MV/cm or 1V/cm as shown in the individual photographs. 


Also, all photographs have a horizontal scale of 5 nsec/cm. 


Fig. 6.10(a) Voltage and current waveforms throughout the 
measuring system 
trace #1 - bias voltage pulse at monitor position A in Fig. 6.5 
trace #2 - bias voltage pulse at the diode holder 
trace #3 - current pulse for avalanching diode taken at 
monitor position B in Fig. 6.5 
trace #4 - detected waveform of the sum of the incident and 
reflected R.F. taken at monitor position C 
Te al CaO ao 
For traces #1 and 2, the 50 ohm impedance of the oscilloscope 
channel was connected to the diode holder in lieu of an avalanching 


Impatt diode. 


For trace #3 the current pulse was taken as a voltage across a 
1 ohm disc resistor in series with an avalanching Impatt diode in the 


diode holding unit. 
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Traces #2, 3 and 4 are for the bias voltage pulse at monitor 


position A, shown in trace #1. 


Fig. 6.10(b) Voltage pulse at monitor position A related to 
the resultant current through the diode measured 


at monitor position B. 


The upper set of three traces are the voltage pulses at 
monitor position A, while the lower set of three traces are the current 


pulses at monitor position B. 


Fig. 6.10(c) Typical measurement photograph 

trace #1 - detected waveform of the sum of the incident and 
reflected R.F. 

trace #2 - detected waveform of the reflected R.F. alone 


trace #3 - the reference voltage pulse at monitor position A 
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a Fig. 6.10 (a) 


Fig. 6.10 (b) 


Fig. 6.10 (c) 
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applied to the Impatt diode was gradually increased until a d.c. current 
of 1 ma was reached. As the reverse bias voltage was increased, the 

R.F. detector was calibrated at several values of reverse bias below 
breakdown. Next, the bias voltage pulses were applied at the same 

levels for which the current response was taken earlier. The oscilloscope 
traces representing the output of the R.F. detector were photographed. 

In the photograph, the magnitude of both the reflected R.F. and the 

sum of the incident and reflected R.F. were recorded, along with the 
waveform at monitor position A. These three traces were recorded 

for each specific level of bias voltage pulse used. A typical photo- 


graph is shown in Fig. 6.10(c). 


The experiments were carried out with two different Impatt 
diodes, whose characteristics are given in Table 6.1]. These character- 
istics were determined from capacitance-voltage measurements taken at 
bias voltages below breakdown and at room temperature. The results of 
the experiments are shown in Figs. 6.11 to 6.23. In these plotted 
results, the measured impedance has been transformed via the equivalent 
circuit of Fig. 6.9(c), to obtain the electronic impedance of the semi- 


conductor diode. Also shown in these figures are the associated 


TABLE 6.1 
CHARACTERISTICS OF H.P. 5082-0437 IMPATT DIODES 


| BREAKDOWN CAPACITANCE AVERAGE DEPLETION 
DIODE} VOLTAGE AT AT BREAKDOWN DOPING LEVEL REGION WIDTH 
1 MA 
16 
#] 101 V 0.34 pf 
#8 1102.2 V O27 ant 
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results from the computer model. These associated results are given 


here for ease of comparison in the following Chapter. 


The measured current response and electronic impedance 
of diode #1 and #8 at room temperature, for various microwave frequencies 
and d.c. bias levels,are given in Fig. 6.11 to 6.23. <A breakdown 


of these figures is as follows: 


- Figs 6.11 to 6.14 - diode #1, 6.0 GHz, 3 and 4 volt 
bias pulses 

- Figs 6.15 to 6.17 - diode #8, 6.0 GHz, 3 and 4 volt 
bias pulses 

- Figs 6.18 to 6.20 - diode #1, 7.8 GHz, 3 and 6 volt 
bias pulses 

~ Figs 6.21 to 6.23 - diode #8, 7.8 GHz, 3 and 6 volt 


bias pulses 


In all the results of the electronic impedance measurements, the plotted 
data displays very little scattering of points and the resulting curves 


are well behaved. 


The series resistance due to contact losses, the resistance 
of the substrate and the undepleted region of the diode can be rep- 
resented by Ro. The value of R. can be taken from the experimental 
results (46)as the real part of the diodes electronic impedance when 
the bias voltage reaches the breakdown value. For the purpose of this 
work the value of No has been taken as the real part of the electronic 


impedance when 1 ma of d.c. current is maintained in the diode. This 
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FIG. 6.11 PULSED D.C. CURRENT RESPONSE FOR IMPATT DIODE #1 AT 6.0 GHz 


oy 50a | 


os | 
Ok Sh 


Chim Bey 

= Al The . 
“he listers tL Baye, be 
a | aaa vi ee | v 


= 


sacar cea, em 2.0 a 
yet wa 7 Bitlet’, | 


145 


EXPERIMENTAL 
COMPUTED 240 


220 


200 


180 


160 


140 


CAPACITIVE REACTANCE (OHMS) 


120 


100 


80 


a LME SI NeNSEC 
FROM F 1G2560 11 


-24 -18 -12 =6 0 +6 vale 
RESISTANCE (OHMS) 


FIG. 6.12 COMPUTED AND MEASURED ELECTRONIC IMPEDANCE FOR DIODE 
#1 AT 6.0 GHz - 4 V EXPERIMENTAL BIAS 
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FIG. 6.14 COMPUTED AND MEASURED ELECTRONIC IMPEDANCE FOR DIODE 
#1 AT 6.0 GHz - 3V EXPERIMENTAL BIAS 
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FIG. 6.16 COMPUTED AND MEASURED ELECTRONIC IMPEDANCE FOR DIODE 
#8 AT 6.0 GHz - 4 V EXPERIMENTAL BIAS 
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FIG. 6.17 COMPUTED AND MEASURED ELECTRONIC IMPEDANCE FOR DIODE 
#8 AT 6.0 GHz - 3V EXPERIMENTAL BIAS 
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FIG. 6.19 COMPUTED AND MEASURED ELECTRONIC IMPEDANCE FOR DIODE #1 
AT 7.8 GHz - 6 V EXPERIMENTAL BIAS 
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value is indicated on all the appropiate figures. For diode #1, the 
value of R. was about 1.5 ohms from impedance measurements taken at 
both 6.0 and 7.8 GHz, while for diode #8 the value of R. was 1.2 ohms 


from measurements at 6.0 GHz,and 2.4 ohms for measurements at 7.8 GHz. 


The value of R. was not included in the computer model 
consequently, when comparing the computed and experimental results, 
the value of R. should be algebraically added to the computed results. 
This is readily accomplished as R. only affects the real part of the 
calculations and the computed curves are simply shifted along the real 


axis an amount equal to the magnitude of R.- 


6.4.3 Accuracy of Measurement System 


Because of the unknown influence of the second harmonic on 
the measurements and the absence of facilities to conduct highly accurate 
steady-state impedance measurements, it is difficult to arrive ata 
specific evaluation of the measurement accuracy for the method used 
here. Further, a detailed investigation of the inaccuracies in the measure- 
ment system would be overshadowed by the limitations inherent in the 
coaxial four-port circulator and R.F. detector used. However, a qualitative 


discussion should prove useful. 


In the Carrier Reinsertion method there are three main areas 
which are subject to measurement errors. These sensitive areas are; 
the packaged diode mounting unit, the separation of the incident and 
reflected R.F. signals and the R.F. detector itself. To provide a 
means of calibrating the R.F. detector sub-system, the impedance of the 


packaged diode should be known for several values of reverse bias voltage 
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below breakdown. The diode mounting unit has to reproduce these diode 
package impedances as magnitude and phase in a reflected R.F. voltage. 
To minimize possible errors being introduced during this transformation, 
the packaged diode holding unit should be an integral part of the 
measurement of the packaged diode impedance, both external to and 
internal to, the Carrier Reinsertion Method. By designing the diode 
holder to be compatible with a standard APC-7 precision connector, this 


source of error was minimized for the work carried out here. 


An important area closely associated with the diode holding unit, 
that may introduce some error in the impedance measurements, is the four- 
port circulator used to separate the incident and reflected R.F. signals. 
In the coaxial four-port circulator used for this work, the isolation 
between the paths followed by the incident and reflected R.F. signals was 
found to be 20 dB. Therefore, it would be possible to have up to 1% of 
the incident R.F. power, or 10% of the incident voltage, bypassing the 
diode holder and interfering with the reflected R.F. signal. As long as 
the magnitude of the reflected R.F. signal is not small with respect to 
the magnitude of the incident R.F. signal, 20 dB of isolation should 


therefore provide experimental results of sufficient accuracy. 


Another possible source of measurement error lies in the operation 
of the R.F. detector sub-system. As discussed in the Appendix B, 
the actual detector diodes have R.F. - input- d.c. - output character- 
istics that vary appreciably when used over a wide range of incident 
R.F. signal powers. This necessitates a piece-wise linear approximation 


to these characteristics and the determination of the constants of 
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proportionality (to calculate the reflection coefficient) for each 
linear segment. To minimize the introduction of measurement error in 
this procedure, the R.F. detector unit should be selected, or tuned, 
for best linearity over the required range of operating levels. Also, 
there should be sufficient versatility in the measurement system to 
ensure that the range of detected output levels evident in the actual 
measurements, was covered by the calibration procedure carried out at 


reverse bias voltages below breakdown. 


The error associated with the Carrier Reinsertion Method of 
impedance measurement should be small and will mainly be determined by 
the characteristics of the R.F. detector unit. Some error was inherent 
in this measurement method, as it was used here, but, this was essentially 
due to availability of equipment and is not believed to have significantly 


affected the overall accuracy of these measurements. 


6.5 Summary 


The measurement of microwave impedance under transient 
conditions requires that the measurement method be capable of following 
variations in the magnitude and phase, of the impedances being measured, 
as a function of time. In this chapter, two distinct types of microwave 
impedance measurement methods have been evaluated, with regard to their 
suitability for the transient impedance measurements on Impatt diodes. 
These measurements of the complex impedance were carried out using a 
direct comparison of the voltage/current waveforms and also, by measuring 


the complex reflection coefficient. 


In the measurement of the transient impedance by the direct 
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comparison of voltage and current waveforms, two methods were used. 
Both the Cavity Oscillator and the Air-line Tee Methods were tested. 
These methods were readily implemented, as a small amount of equipment 
was required, however, both these measurement methods proved to have 


serious limitations that restricted their usefulness. 


The Carrier Reinsertion Method for the measurement of transient 
impedance actually measures the complex reflection coefficient, which 
is then converted to an impedance. This method requires more extensive 
use of test equipment but is free of the difficulties associated with 
the waveform comparison methods. To the best of this author's know- 
ledge, the Carrier Reinsertion Method, as developed here, has not been 
used in any previous impedance measurements on Impatt diodes. This 
method is suitable for either steady-state or transient impedance measure- 
ments. The method, as utilized here, had some limitations due to 
availability of equipment but overall, the method appears to have good 


potential for future work in this field. 


The electronic impedance of Impatt diodes, under transient 
conditions, has been measured by using the Carrier Reinsertion Method. 
~The plotted data form well behaved curves that display very little 
scattering of the measured values. These measurements have been obtained 
by reverse biasing the Impatt diode to maintain a d.c. current of 
1 ma and then modulating the d.c. Mine voltage supply with low-voltage 
pulses. These bias pulses had rise times as short as 5 nanoseconds 
at the reference plane of the diode and it was determined that the d.c. 


current response of the packaged diode was at least equivalent to these 
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bias voltage pulse rise times. Therefore, short duration bias pulses, 
in the range of 30 nanoseconds, could be used. These short duration 
bias pulses enabled the use of very low duty cycles which minimized 
the effects of temperature changes in the diode's semiconductor junction. 
The measurement holding unit for the packaged diode was designed to 
give a simple transformation from the measured impedance values to 

the electronic impedance of the semiconductor chip itself, thus, 
simplifying the calculation and improving the accuracy of the measure- 
ments. Since these impedance measurements have been referenced to the 
d.c. current response of the diode, comparisons of these measured 
results with computer models of the Impatt diode are readily carried out. 


A detailed comparison will be given in the following chapter. 
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CHAPTER VII 
Discussion of Theoretical and Experimental Results 


7.1 Introduction 


US3et e147 040) Geer bind 


There have been several publications 
theoretical models of the Impatt diode, but very few publications 
include a direct comparison between the theoretically-derived and the 
experimentally-measured electronic impedance. The work undertaken 


here has been carried out with the intent of substantiating the 


validity of the theoretical transient model with experimental results. 


In the experiments, and hence in the theoretical model, 
the bias voltage pulses were of short duration (about 30 nanoseconds). 
For these short-durations bias pulses, temperature effects on Impatt 
diode operation are minimal. Further, the comparison between the 
computations and the measurements were made on the basis of similar 
d.c. current pulses and thus, changes in the temperature of the semi- 
conductor junction are accounted for. Therefore, the discussion of 
the transient impedance characteristics will be carried on without 


any reference to the junction temperature. 


7.2 Comparison of Computations and Measurements 


The experiments were carried out on two separate diodes and 
at two frequencies that were sufficiently distinct that differences 
in the behaviour of the electronic impedance (20) could be observed. 
This procedure ensured a more valid comparison between the computed and 


the experimental results. As an indication of the type of comparison 
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to be expected, it was noted that Greiling and Haddad! 2) in their 
comparison between theory and experiment, encountered differences in 
the real and imaginary part of Zo in the range of 15-60% and 5-102, 
49) 


respectively. Further, Decker ée¢ au! claim that an accurate 
comparison between theory and measurements can only be obtained when 


the effects of injecting contacts are taken into account. 


Generally, the computations and the measurements of the 
diode's electronic impedance exhibited the same type of behaviour. 
This agreement holds for both frequencies used and also at high and 
low bias levels. A comparison of the plotted results in Figs. 6.11 to 
6.23 shows that the slopes and curvatures of the computed and ex- 
perimental results are very similar. A significant indicator for 
the comparison of the results is in the behaviour of the electronic 
impedance when the bias pulse passes from its initial rapid rise to 
its relatively flat top. Here, the clockwise loop at 6.0 GHz, Figs. 
6.14 and 6.17, and the counter clockwise loop at 7.8 GHz, Figs. 6.19 
and 6.23, are displayed by both the computed and the experimental 
results. About the only significant discrepancy between the computed 
and the experimental results appears to be in the relative magnitudes 
of the real and imaginary parts of the electronic impedance, with the 
largest discrepancy occurring in the real part. The differences between 


the computations and the measurements are typically 20% for diode #8 


and 40% for diode #1 in the real parts of Zoe with the differences in the 
imaginary parts of Ze being typically 8% for both diodes. When making these 
comparisons, it should be remembered that the plotted computed results have 


not taken the value of R. into account. (Where R. is the resistance 
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due to ohmic contacts, substrate resistance and undepleted portion of the 


drift region. Re is in series with Za): 


From a comparison of the plotted results, it is evident that 
the measurements, in some aspects, do not exhibit the same magnitude of 
variation in behaviouras do the computed results. This is probably due 
to a slower response in the measuring system, although part of the 
magnitude discrepance could be due to the influence of the second 
harmonic in the driving R.F. voltage waveform, as it appears across the 


diode. 


7.2.] Response of the Measuring System 


From Figs. 6.12, 6.19, and 6.22, the experimental results 
appear to be a smoothed version of the theory. Part of this difference 
between the computations and the measurements could be due to the response 
in the R.F. detection subsystem. This detector response could account 
for the discrepancy between the computed and measured behaviour of the 
electronic impedance during the decay time of the bias voltage pulse. 
During the initial part of the decay time, both the computed and measured 
results show the same relative shift in the values of Le but, as the 
decay time progresses the computations show a recovery trend that is 
not apparent in the experimental results, until after a few nanoseconds 
have passed. 

Equipment problems alone, do not account for all the differences 


between the computation and the measurements. This is evident in that 


the two diodes are affected differently, with the measured results for 


77 oA oe. 
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diode #1 displaying a more restricted behaviour than those of diode #8. 
This is particularly evident at the higher current levels, PTS 0 tLe 
and 6.16. At the lower bias levels the behavior of the two diodes is 


much closer related. 
7.2.2 Second Harmonic Effects 


The current waveform of the diode computer model is Fourier- 
analyzed to obtain the d.c. component plus, the first, second and third 
harmonics. The ratio of the magnitude of the second harmonic to the 
magnitude of the fundamental frequency component (Yo/1) is plotted in 
Figs. 7.1 and 7.2 for 6.0 and 7.8 GHz, respectively. Generally, there 
iS an increase in the value of Yoyq? with increasing d.c. current level 
in the diode and the overall behaviour of the Y2/] plot is similar 


to the waveform of the corresponding d.c. current pulse. 


For diode #8, the value of Yo/] is relatively small and varies 
from 1 to 4%, with different frequencies and bias voltage levels. However, 
for diode #1, the value of Yo/] shows a Significant change for different 
bias voltage levels, ranging as high as 14% for high d.c. current levels 
at 6.0 GHz. For both diodes the value of Yo/1 is larger for 6.0 GHz 


than for 7.8 GHz with the biggest difference being exhibited for diode 


#1. 
(25) 
Greiling and Haddad have shown that it is possible for an 


appreciable second harmonic component to be present in the R.F. voltage 
waveform, generated by the second harmonic component in the R.F. current 
waveform of the diode. In the experiments carried out here, the driving 


R.F. voltage applied to the Impatt diode is initially of single frequency. 
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FIG. 7.1 COMPUTED RATIO OF SECOND HARMONIC TO FUNDAMENTAL R.F. 
COMPONENT OF IMPATT DIODE CURRENT AT 6.0 GHz 
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FIG. 7.2 COMPUTED RATIO OF SECOND HARMONIC TO FUNDAMENTAL R.F. 
COMPONENT OF IMPATT DIODE CURRENT AT 7.8 GHz 
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Then, since the second harmonic component of the R.F. voltage is not 

absorbed in the measuring system (coaxial four-port hybrid, Fig. 6.1), 
it is possible that a significant R.F. voltage at the second harmonic 
frequency could be developed in the diode holder, as the Impatt diode 
is biased to higher current levels. Consequently, another variable is 


introduced to modify the measured values of the electronic impedance. 


A possible explanation of the significant differences in the 
magnitude of the electronic impedances, of diode #1, with respect to the 
computed results, at the higher current levels, could be the influence 
of a second harmonic component in the R.F. voltage. From Figs. 6.14 
and 6.20, it is observed that the measured and computed values of the 
electronic impedance, for diode #1, display a close correspondence for 
both 6.0 and 7.8 GHz. These results are for moderate current levels in 
the diode. For higher currents the results are as shown in Figs. 6.12 
and 6.19, where there is significant differences between the computed 
and measured values of the electronic impedance. From Fig. 7.1 and 
7.2, the value of Yo/] is much larger for the higher current levels. 
This discrepancy between the computed and measured values of the elect- 
ronic impedance, at the higher current levels, is not evident for 
diode #8. Also, the value of Yo? is relatively small for both moderate 


(50) shat the 


and high current levels in this diode. As it is well known 
second harmonic in the R.F. voltage does have a significant influence 
on the behaviour of the Impatt diode, it is likely that the second harmonic 


is also influencing the measurements taken here . 
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7.2.3 Degree of Correspondence of Computed and Experimental Results 


Because of the questions raised concerning the response of 
the measuring system and the extent of the effects of the second 
harmonic component of the R.F. voltage across the diode, it was thought 
that only a limited amount of curve fitting should be carried out. 
Therefore, efforts to fit the computations to the measurements were 
confined to small variations in the values of the parameters Xo Ke Cy 
and Tq in the computer model. The effects of variations in these 
parameters are outlined in Chapter IV, sub-section 4.3.2, as are the 
final values chosen for the comparison between the experimental and 


computed results (Table 4.1). Consequently, only the significance of 


the changes from the theoretical values will be discussed here. 


The value of the equivalent width of the avalanche region, 
X, > was arrived at by considering the experimental and computed 
characteristics of diode #1. From Fig. 4.6, it is evident that the value 
of x, is less than 155 microns while, from the characterization of 
Scharfetters model (Fig. 3.4), the value of X 4 is at least 1.35 microns. 
Since the doping levels oF diode #1 were slightly less than those 
of Scharfetters model, a value of 1.45 microns for X,° of diode #1, was 
selected. Diode #8 has a doping concentration that is less than that 
of diode #1 and therefore has a value of Xs greater than 1.45 microns. 
From information <> given on the predicted widths of avalanche regions 
for different doping levels of Silicon, a value of 1.65 microns for Xa? 


in diode #8 was deemed to be appropiate. The values of Xs chosen for 


these diodes are in general agreement with those predicted by previously 
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published work and were in proportion to the values of K selected for the 


diodes. 


It is difficult to correlate the sum of the avalanche and 
drift regions with the total depletion width of the diode. The value 
of x, is an equivalent width, while the width of the drift region and its 
associated transit time are actual values. Therefore, the value of Ty 


was chosen as outlined in Chapter IV and then held constant. 


Magnitude scaling was mainly accomplished by variations, in 
the values of the depletion layer capacitance (C,). The variations from 
the measured values of Cy at 1 MHz were a maximum of 0.03 pf and this 
is within the tolerance between capacitance values measured at 1 MHz 
and several GHz, when stray capacitances of the package are taken into 
teat tana Therefore, the values of Cy used in the final comparison 


between computed and experimental results are quite accurate and within 


the tolerances found by other workers in this field. 


As outlined in Chapter IV, the value of the coefficient for 
the space charge term associated with the width of the avalanche region 
(K), was selected to provide the best behavioural comparison between 
the computations and the measurements. Some magnitude scaling was also 
achieved by this selection but, this was not the main effect. From 
the constant field theory, the values of K should be 0.011 for diode #1] 
and 0.015 for diode #8. As shown in Table 4.1, the actual values of 
K used differed a little from these theoretical values and also displayed 
some frequency dependence. Some differences between the theoretical 


and actual values used are to be expected, as the constant field model 
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only gives an estimate of the true value of K. The different values of 

K for different frequencies would seem to indicate a functional 
dependence for K that is not evident in the constant field model. 
However, this frequency dependence is small and may have been accentuated 


by other characteristics of the measurement system. 


7.3 Effects of the Bias Pulse on the Electronic Impedance. 


As described in Chapter IV, a comparison between computed 
and experimental results must be made on the basis of the similar 
d.c. current pulses, which result from the simulated and applied bias 
voltage pulses, respectively. Both the computations and the measure- 
ments indicate that this resultant d.c. current pulse will be very 
Similar to the bias voltage pulse, across the diode, in rise times and 
overall shape. As the general shape and rise time of the bias voltage 
pulse, across the diode, can be controlled by the circuit designer, 
advantage can be taken of the dependence of the electronic impedance 


upon the d.c. current pulse in the diode. 


7.3.1 Magnitude of the Bias Voltage Pulse 


Generally, as the magnitude of the d.c. current increases, 
the magnitude of the negative impedance of the diode increases as well. 
This increase contributes to an increase in the R.F. power output of the 
diode. However, from Figs. 6.12 and 6.16 it is observed that, once the 
level of d.c. current exceeds a certain value, the real part of the 
electronic impedance starts to decrease in magnitude. This is a result 
of the space charge effect of the mobile carriers modifying the total 


phase shift across the diode. The larger the space charge effect, the 
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shorter the time allowed for the avalanche current pulse to form. This 
results in a smaller phase shift being developed in the avalanche region, 
between the conduction current and the R.F. voltage. The change 

from increasing magnitude to decreasing magnitude, in the real part of 
the electronic impedance, is clearly evident at 6.0 GHz but, is not 
evident at 7.8 GHz, Figs. 6.19 and 6.22, even though the current levels 
at 7.8 GHz were much higher than those at 6.0 GHz, Figs. 6.11 and 6.18, 


respectively. 


The reason for the reduced effects of the space charge at 
7.8 GHz is that, at this higher frequency, the phase shift across the 
avalanche region is a smaller part of the total phase shift across the 
diode. Taking a mean length for the drift region as 3 microns for the 
diodes used here, at 6.0 GHz this length represents about 65° of 
phase shift, between conduction current and R.F. voltage, while at 7.8 
GHz this same length represents about 86° of phase shift. Consequently, 
the reduction of the phase shift in the avalanche region, due to the 


Space charge effect, is much less critical at 7.8 GHz than at 6.0 GHz. 


7.3.2 Rise Times of the Bias Voltage Pulse 


The effect upon the electronic impedance of the rise of the 
bias pulse can be obtained from the plotted results of the electronic 


impedance as a function of time and d.c. current levels. The plotted 


— 


results were analyzed and salient features are summarized in Tables 7. 
and 7.2; the information was arrived at in the following manner: 
consider diode #1 at 6.0 GHz; selecting a current level of 30 ma in 
Figs. 6.11 and 6.13, the corresponding times are 3 nanoseconds and 


4 nanoseconds in Figs. 6.12 and 6.14 respectively. The values of the 
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TABEE"7 <1 


TRANSIENT AND STEADY-STATE VALUES OF ELECTRONIC IMPEDANCE 
COMPUTED AND MEASURED AT 6.0 GHz 


DIODE | FREQ |METHOD j|CURRENT 3 VOLT 6 VOLT STEADY 
(GHz) LEVEL BIAS BIAS STATE 


Computed |Fig.6.11 | Fig.6.14 Fig.6.12 Fig.7.4 
4 nsec 3 nsec 
-13.6-j106 -14.5-j107.5 | -12.8-j104 
: 12 nsec 
-18.3-j144.5 -18.5-j146 
4.7 nsec 
-20.8-j148 -18.5-j147 
7 nsec 
-6.7-j205 -4,0-j205 
12 nsec 
$y ier ao +6.4-j225 
Measured s 3 nsec 
4-j -7.5-j99 


Pag’ Om iy: Fig. 6. 16 

4 nsec 3 nsec 

-12.6-j113 -12.8-j113.5 | -11.6-jl12 

12 nsec 4.8 nsec 

-15.4-j135 -17.2-j139 -15.4-j135 
7 nsec 
-13.2-j159 -12.5-j158 
15 nsec 
-8.5-j169 -8,4-j169.5 

Measured|25 ma 3.75 nsec 3.5 nsec 
-10.4-j120 -10.8-j124 
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TABLE 7.2 
TRANSIENT AND STEADY-STATE VALUES OF ELECTRONIC IMPEDANCE 
COMPUTED AND MEASURED AT 7.8 GHz 


DIODE}FREQ | METHOD CURRENT | 3 VOLT 6 VOLT STEADY 
(GHz) LEVEL BIAS BIAS STATE 


#1 17.8 Computed Fid.6.20 Pig.oe19 Fide as 
4 nsec 2.3 nsec 
-9.3-j64 -7.0-j67.3 -8.2-j68.0 
6 nsec 
-12.4-j76 -12.6-j75 
23 nsec 
-12.6-j75 
10 nsec 
-14.5-j78 -14.4-j78.5 
6 nsec 
-27.8-j113 -28-jll2 
12 nsec : 
-3].6-j134 -32.2-j131 
17 nsec 
-32.6-j135 
15 nsec 
-33.8-j142 -32.8-j137 
Measured 2.5 nsec 
-3.5-j74.5 
| 


Computed Eig). 6823 Fidg..6%22 
5 nsec 3 nsec 
-11.4-j85 -10.6-j83 -10.2-j84.5 
9 nsec 


-13.0=387 8 -12.9-j87.5 | 


6 nsec 

-22.5-j108 -22.4-j106 

10 nsec 

-25.4-j113.5 -25.0-j113.5 

15 nsec 

-26.0-j117.0 -26.0-j117 
Measured 3.25 nsec 
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electronic impedance are then taken from the appropiate graphs. Values 
of the computed steady-state electronic impedance are taken from 

Figs. 7.3, 7.4 and 7.5 for the appropiate d.c. current levels. These 
Steady state values of electronic impedance were calculated using 

the same parameter values as those used for the final computations 


plotted in Figs. 6.12 to 6.23. Thus, a direct comparison is possible. 


Effects of the mobile carrier space charge will place an 
upper limit on the magnitude of the real part of the electronic impedance 
(Z,). Apart from this, the rise times of the bias voltage pulse can 
have a significant effect upon the characteristics of the diodes Zee 
Essentially, the effects of the rise times are confined to modifying 
the real part of Ze with very little change in the imaginary part of 


Le 
e 


From the data given in Tables 7.1 and 7.2, it is evident 
that the maximum negative value of the real part of Ze can be increased 
beyond the equivalent (same d.c. current level) steady-state values 
for both diodes and at both frequencies used here. However, there is 
a limit beyond which a further decrease in the rise time of the bias 
voltage pulse produces little further increase in the negative real part 
of Ls and may even reduce the magnitude of the real part of Za This 
effect is well illustrated in Table 7.2 for 7.8 GHz and for 6.0 GHz in 
Table 7.1 and Figs. 7.6 and 7.7, with the corresponding current levels 
being given in Figs. 6.13 and 6.15, respectively. This minimum rise time 
appears to be about 5 nanoseconds for 7.8 GHz and 4 nanoseconds for 


6.0 GHz. It would appear that by using the minimum rise time of the 
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FIG. 7.3 COMPUTED STEADY-STATE VALUES OF ELECTRONIC IMPEDANCE 
AT 7.8 GHz AND 2 VOLTS R.F. 
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FIG. 7.4 COMPUTED STEADY-STATE VALUES OF ELECTRONIC 
IMPEDANCE AT 6.0 GHz AND 2 VOLTS R.F. 
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FIG. 7.5 COMPUTED STEADY-STATE VALUES OF ELECTRONIC IMPEDANCE 
AT 6.0 GHz AND 10 VOLTS R.F. 
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FIG. 7.6 COMPUTED VALUES OF ELECTRONIC IMPEDANCE FOR DIODE #1 
AT 6.0 GHz FOR DIFFERENT BIAS VOLTAGE PULSE RISETIMES 
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FIG. 7.7 COMPUTED VALUES OF ELECTRONIC IMPEDANCE FOR DIODE #8 
AT 6.0 GHz FOR DIFFERENT BIAS VOLTAGE PULSE RISE TIMES 
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bias voltage pulse an increase in the magnitude of the real part of 
Le of about 10-15% can be expected during the rise time of the bias 
voltage pulse. This increase is not very large but may be useful in 


the start of oscillations in a coupled cavity. 


7.4 Recovery Time from Transient to Steady-State Conditions 


Overall, there are relatively small departures from the 
steady-state values of the electronic impedance during transient 
conditions. These departures occur mainly during the buildup and decay 
of the bias voltage pulse, with the transient behavior of Ze during the 


decay, being more pronounced and of longer duration. 


From the appropiate Figures and Tables 7.1 and 7.2, it is 
apparent that the effects of the rise times of the bias voltage pulse 
are somewhat different at 6.0 GHz from those at 7.8 GHz. This difference 
is probably due to the fact that the real part of L., reaches its 
maximum magnitude at relatively low current levels at 6.0 GHz, while 
at 7.8 GHz this maximum limit has not quite been reached at very high 
current levels. The clockwise loops at 6.0 GHz, (Figs. 6.14 and 6.17), 
and the counter-clockwise loop at 7.8 GHz, (Figs. 6.19 and 6.23), in 
the behavior of Zee are due to the transient values of L. returning to 
their steady-state values. At 6.0 GHz, during the transition interval 
of the d.c. current pulse from its rapid rise to its relatively slowly 
varying top, the more negative values of the real part of the transient 
Z, monotonically approach their steady- state values; hence, the clock- 
wise loop. For 7.8 GHz, during the same transition of the d.c. current 


pulse, the transient Le appears to overcorrect, with the real part 
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FIG. 7.8 COMPUTED CURRENT RESPONSE FOR DIODE #1 AT 6.0 GHz 
WITH 10 VOLTS R.F. AND 4 VOLTS EXPERIMENTAL BIAS 
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FIG. 7.9 COMPUTED VALUES OF ELECTRONIC IMPEDANCE FOR DIODE #1 
AT 6.0 GHz AND 4 V EXPERIMENTAL BIAS 
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FIG. 7.10 COMPUTED CURRENT RESPONSE FOR DIODE #1 AT 6.0 GHz 
WITH 10 VOLTS R.F. AND 3 V EXPERIMENTAL BIAS 
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FIG. 7.11 COMPUTED VALUES OF ELECTRONIC IMPEDANCE FOR DIODE #1 
AT 6.0 GHz AND 3 VOLTS EXPERIMENTAL BIAS 
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FIG. 7.12 COMPUTED CURRENT RESPONSE FOR DIODE #8 AT 6.0 GHz 


EXPERIMENTAL 


COMPUTED 10 V R.F. 


6 12 18 PACE FOR 030 
TIME (NSEC) 


WITH 10 VOLTS R.F. AND 3 VOLTS EXPERIMENTAL BIAS 


186 


36 


JATMAMTRIAN 
AA V OF aaTuaMOD 


187 


% COMPUTED 2 V R.F. 
2 Se GOMPUTED 10 VERE. 


140 


130 


120 


110 


100 


Oy ys ae TLME GENE NSEC FROM 90 
FIGS.7.12 AND 6.15 


-20 -16 -12 -8 -4 0 +4 
RESISTANCE (OHMS) 


FIG. 7.13 COMPUTED VALUES OF ELECTRONIC IMPEDANCE FOR DIODE #8 
AT 6.0 GHz AND 3 VOLTS EXPERIMENTAL BIAS 
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FIG. 7.14 COMPUTED CURRENT RESPONSE FOR DIODE #8 AT 6.0 GHz 
WITH 10 VOLTS R.F. AND 4 VOLTS EXPERIMENTAL BIAS 
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FIG. 7.15 COMPUTED VALUES OF ELECTRONIC IMPEDANCE FOR DIODE #8 
AT 6.0 GHz AND 4 VOLTS EXPERIMENTAL BIAS 
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becoming less negative than its steady-state value and hence the 
counter-clockwise loop results. In both cases, the time required for 
the transient Le to settle to their steady-state values is about 5 


nanoseconds. 


The previous discussions have been for the case where the 
magnitude of the R.F. voltage has been one-or-two volts across the 
diode. To gain an insight into the behaviour of the electronic impedance 
at higher R.F. voltages, some computations were carried out using 10 
volts as the magnitude of the driving R.F. voltage across the diode. 
These computations were restricted to 6.0 GHz as the bias levels used 
at this frequency give practically all the behavioural aspects of the 
electronic impedance for the diodes under study here. The results of 
these computations, for diode #1 and #8 are given in Figs. 7.8 to 7.15 
inclusive. As expected, the magnitudes of Z. calculated using ten 
volts R.F. are less than those using two volts R.F. However, a 
comparison of the transient and steady-state values of Zee under the 
influence of 10 volts R.F., shows very little difference between the 
two sets of computations. Therefore, a discussion of the recovery time 
from transient to steady-state conditions is somewhat unimportant and 
for practical purposes, it can be assumed that steady-state conditions 


prevail at all times for this higher magnitude of R.F. voltage. 


7.5 Summary 
Generally, there was found to be a close comparison in the 
behavioural aspects of the computed and measured values of the electronic 


impedance under transient conditions. Comparisons between the computations 
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and the measurements were carried out on the basis of similar d.c. 
current responses for the simulated and actual bias voltage pulses 
used. This, plus the short-duration bias pulses, ensured that the 
effects of changes in the temperature of the semiconductor junction 


could be ignored. 


The main discrepancy between the computed and the measured 
results appears in the relative magnitude of Les An important cause 
of this discrepancy is due to the influence of the second harmonic 
component of the R.F. voltage impressed upon the diode. The ratio 
of the magnitude of the second harmonic to the fundamental in the R.F. 
current waveform of the diode, increases with increasing d.c. current 
levels. A second harmonic component in the R.F. voltage is developed 
from this second harmonic in the R.F. current and as the measuring 
system will not absorb this second harmonic frequency, it is likely 
that a significant second harmonic component is present in the R.F. 
voltage in the diode holder. A detailed evaluation of this second 


harmonic effect was not carried out in this study. 


In general, the significant differences between transient 
and steady-state values of the electronic impedance, occur during the 
initial buildup and decay of the bias voltage pulse across the diode. 
These differences are most prominent at the lower magnitudes of R.F. 
voltage. At R.F. voltages of one or two volts, the transient values 
of the electronic impedance can have real parts that are about 10-15% 
more negative than their steady-state counterparts. Under these 


conditions, there appears to be an optimum rise time of 4 to 5 nano- 
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seconds for the bias voltage pulse, with a recovery time of about 5 
nanoseconds needed for the transient values of the electronic 
impedance to return to their steady-state values. At magnitudes of 
R.F. voltage in the range of 10 volts, there is negligible difference 


between the transient and steady-state values of the electronic impedance. 
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CHAPTER VIII 
DISCUSSION AND CONCLUSIONS 


This study has investigated the behaviour of the electronic 
impedance under transient conditions for a type of Impatt diode. The 
analysis and results are applicable to Read diodes and one-sided 
abrupt p-n junction diodes. In the course of this work, several 
aspects of the operation of these diodes have been investigated, with 


the following results: 


Th: An improved analytical model of the diode, which 
applies to both steady-state and transient operating 
conditions, has been formulated and evaluated. 

rr A model for the transient heat flow in the diode's 
semiconductor junction has been developed which has 
been verified by experimental results and incorporated 
into a computer model of the diode. 

os A method for the measurement of the diode micro- 
wave impedance has been developed which allows the 
measurements to be conducted under steady-state or 


transient conditions. 


All of these aspects have been incorporated into a computer model of the 


Impatt diode which has been programmed on a large digital computer. 


The analytical equations developed here describe the behaviour 
of the Impatt diode under d.c., small-signal and large-signal conditions 


for both steady-state and transient operation. These equations are 
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relatively simple in form and make use of a minimum number of approx- 
imations, while still maintaining an efficient computer solution. The 
equations are sufficiently versatile to describe either pn or np 
one-sided abrupt junctions in a variety of semiconductor materials. 
The structure of the equations is such that the effects of the several 
pertinent parameters are readily apparent, thus facilitating the 
Simulation of the experimentally obtained results with the computer 
model. This leads to a more rapid understanding of the problem under 


study. 


The describing equations and the computer model for the Impatt 
diode have been shown to be valid from comparisons with other published 
works and from experimental results obtained for the transient behaviour 
of the electronic impedance of the diode. These equations have been 
used to accurately model the steady-state and transient thermal response 
of the p-n junction along with the d.c. and a.c. operation of the diode 
under steady-state and transient operating conditions. An empirical 
representation for the current profiles in the avalanche region of the 
diode,under dynamic operating conditions,has allowed the equations to 
include the individual ionization rates and scattering limited drift 
velocities for the carriers and also has allowed the inclusion of the 


space charge term for the avalanche region. 


The effects of injecting contacts and the spreading of the 
avalanche and drift regions under conditions of high R.F. signal levels 
have not been taken into account. Similarly, the associated effects 


of non-scattering-limited drift velocities, unswept epitaxial layer 
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and dispersion of the pulse of carriers in the drift region have been 
omitted from this analysis. For the model of the Impatt diode under 
Study here, all of these factors are of somewhat secondary importance, 
except for the effects of injecting contacts, which may influence 

the degree of correspondence between the magnitude of the computed 


and measured values of the diode electronic impedance. 


In studies concerning the transient impedance of Impatt 
diodes for applications in oscillators or amplifiers, the behaviour of 
the junction temperature as a function of the duration of the bias 
voltage pulse is important. Both the steady-state and transient thermal 
properties of a type of Impatt diode were investigated. A study of 
the time dependence of the junction temperature was carried out over 
an ambient temperature range of 20°C to 150°C and over a time interval 
of from 5 to 200 nanoseconds after the onset of the bias voltage pulse. 
These studies have related the experimental results to a theoretical 
model for heat flow in a semiconductor junction. Good agreement has 


been obtained between the two types of results. 


A theoretical model of the time dependence of the junction 
temperature was arrived at by considering a one-dimensional analysis 
of a semi-infinite semiconductor with a region of heat generation at 
one end. This model is sufficient for pulsed operation of the diode 
with pulse widths and repetition rates such that the high thermal 
resistance of the bonding wire does not influence the heat flow in the 
substrate. Incorporating the time dependence of the junction temperature 
into the computer model, for the Impatt diode, was readily carried out 


by adding a time varying term to the simulated d.c. bias voltage pulse. 
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This computer model produced very good agreement with the experimental 


results. 


It has been found that the relative change in junction temper- 
ature during the application of the bias voltage pulse, was essentially 
the same for the entire range of heat sink temperatures used here. 

Thus, the thermal behaviourof the Impatt diode will remain essentially 
the same for average junction temperatures from 20°C to 150°C. 
Consequently, the junction temperature does not have to closely approach 
the heat sink temperature between each bias pulse, as only the average 
junction temperature is important and higher pulse repetition rates 


can be used without altering the transient thermal response of the diode. 


The junction temperature mainly alters the d.c. current 
level in the diode, therefore, as long as the computed and experimental 
results are compared on the basis of: similar d.c. currents in the 
diode, the effects of junction temperature changes on this comparison 
are minimized. If the bias voltage Bin ses are of short duration (around 
40 nanoseconds), the effects of junction temperature changes are 
significantly reduced and may even be able to be ignored in certain 


circumstances. 


A method has been developed for impedance measurements under 
transient conditions. In this method, the reinsertion of the R.F. 
Carrier allows the measurement of the complex reflection coefficient 
as a function of time. This Carrier Reinsertion Method is suitable 


for either steady-state or transient impedance measurements. Under 
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transient conditions the response time of the measuring system was 

found to be as low as 5 nanoseconds utilizing standard, commercially 
available,.coaxial and waveguide components. This response time 

can be further reduced by modifying the manner in which the bias 

voltage pulses are applied to the packaged diode. The accuracy of 

this measuring system should be at least equivalent to that of 
conventional VSWR techniques used for steady-state impedance measurements , 
although a detailed investigation of the relative accuracy was not 
carried out. Overall, the Carrier Reinsertion Method appears to have 
good potential in the measurement of microwave impedances that are 


a function of time. 


It has been found that the significant differences between 
transient and steady-state values of the Impatt diode electronic 
impedance (Z,); are of relatively short duration and are confined to 
the initial buildup and decay of the bias voltage pulses applied to 
the diode. These differences are most prominent at the lower R.F. 
voltages and are largely confined to the real part of Zee with the 
imaginary part of rie being relatively unaffected. At R.F. voltages 
of one or two volts, the transient values of the real part of Ze can be 
about 15% more negative than their steady-state counterparts. Under 
these conditions, there appeared to be an optimum rise time of 4 to 5 
nanoseconds for the bias voltage pulse across the diode, with a recovery 
time of about 5 nanoseconds needed for the transient values of Le to 
return to their steady-state values. At the higher R.F. voltages, such 


as 10 volts, there appeared to be a negligible difference between the 
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transient and steady-state values of Zee 


These variations in the real part of Z. are of significance 
when the Impatt diode is employed in R.F. oscillators. By utilizing 
fast rise time bias voltage pulses, the initial growth of oscillations 
in a coupled cavity will be enhanced by the more negative real part 
of Zee Further, it appears that this enhancement is readily attainable, 
as rise times of the bias voltage pulse across the diode of about 
4 to 5 nanoseconds, can be obtained with commonly available microwave 
circuit components. The relatively short duration of the transient 
effects ensures that the steady-state values of Le can be used for most 
design purposes. In self-deviating Impatt oscillators, used in 
frequency modulated carrier systems, the information is placed on the 
carrier by modulating the bias supply of the Impatt diode. When pulse 
code modulation is employed some additional amplitude distortion may 
be introduced due to the transient behaviour of the real part of Ze: 
However, the magnitude of the frequency deviation will not be affected 


by the transient conditions. 


An interesting result of this investigation of transient 
conditions on the electronic impedance of Impatt diodes is that the 
real part of L may become positive under certain conditions of bias 
voltage and diode design. It appears, for Impatt diodes which have 
a drift region transit time that is less then 1/5 the period of the 
R.F. frequency of operation, that the real part of Z will become positive 
for moderate, pulsed, d.c. current levels. Under these conditions, 


the maximum negative value of the real part of Le is reached for 
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relatively low d.c. current levels for pulsed operation. When the 
R.F. frequency is increased such that the drift region transit time 
is about 1/4 of the R.F. period, the maximum negative value of the 
real part of La is not reached until relatively high levels of pulsed 


d.c. current are applied. 


This "peaking" characteristic of the real part of Za: with 
increasing d.c. current levels, may be used to advantage in self- 
deviating Impatt oscillators. When the bias circuit of the Impatt is 
modulated, the R.F. output of the oscillator varies in frequency and 
magnitude due to the changes in the imaginary and real parts,respectively, 
of the diode Zee The magnitude variations in the R.F. are undesirable 
and can be minimized by biasing the Impatt diode such that the average 
value of the real part of Ze is at its peak or maximum negative value. 
This minimizes the variations in the real part of Le and maximizes 


the variations in the imaginary part of Za 


This peaking of the real part of the Impatt diode's Ze is 
detrimental to the operation of R.F. amplifiers and pulsed R.F. 
oscillators. As the magnitude of the R.F. voltage across the diode 
increases, less d.c. bias voltage is required for a given level of d.c. 
current. In large-signal amplifiers, the high magnitude of R.F. signals 
may drive the real part of the diode's ZA past its maximum negative 
value, which will result in less gain and may further complicate the 
associated compensation circuitry. For Impatt oscillators that provide 
pulsed R.F., this peaking effect of the real part of Ze will place 


limitations on the permissible magnitudes of d.c. bias voltage and the 
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resulting R.F. signal levels. 


The overall effect of this peaking of the real part of the 
Impatt diode's electronic impedance will be to place more stringent 
limitations on the design of the diode for its specific task, as it 
is the width of the drift region which determines the lower frequency 


and level of d.c. bias current for which the peaking takes place. 


During the course of this work, certain aspects of the 
study arose which merit further investigation: One important topic is the 
influence of the second harmonic of the R.F. voltage across the diode. 
This second harmonic effect limits the validity of the experimental 
results to those obtained with low values of the fundamental R.F. 


voltage. 


Another topic that warrants further investigation is the 
effect of injecting contacts. It has been shown that injecting 
contacts can influence the comparison of computed and measured results 
of the diode's electronic impedance but, the extent of this influence 


has not been established. 


The results of the suggested investigations, combined with 
those of previously published work and the results reported here, should 
enable an analytical model of the Impatt diode to be constructed which 
will produce much closer correspondence between calculated and measured 


values of the electronic impedance. 
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APPENDIX A 
FOURIER ANALYSIS AND COMPLEX IMPEDANCE 


A series of the form 


A, [A, cos(nx) + B, sin(nx)] (1) 


is called a trigonometric series. It is called a Fourier series if 


all the coefficients Ane B, can be obtained in the following manner 


oT 
A, =2 i f(x) cos(nx) dx n=0,152,..- (2) 
O 
oT 
B. -1 f f(xesin(nx) dx. . nel42,34,.- (3) 
O 


By using a change of variables 
T 
=i t 
soa 
the series is now written in a more convenient form as 


f(t) = A/2 + j 


a Ok aes (CE eee EL sthh (Caesoh a (4) 


L 


where the coefficients are given by 
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and 2L is the period over which the integration takes place. 


From Section 4.2 Eq. (18), the total current in the 


external circuit of the Impatt diode can be written as 


Tint see) 2, alee (7) 


d oT 
The external conduction current is a well behaved function and both 
integrable and differentiable over any one period. Therefore, it may 
be assumed that this function has a Fourier series. Similarly, the 


R.F. voltage may be assumed to have a Fourier series also. 


To arrive at the electronic impedance of the semiconductor 
diode, the complex ratio of the fundamental components of Eq. (7) 
and V(t) must be taken. The fundamental component of current may be 


obtained by Fourier analysis of Eq. (7) and written as 
2 a ny fee 
Ir (t)¢ = A cos ( L t JG Bas inet tm) (8) 


This equation is now rewritten in a more convenient form by using 


Euler's formula 


exp ( j rt ) = cos( me ) + j sin ( ci ) (9) 
After some rewriting, Eq. (8) becomes 

Ty(t)p = Al + Bo sin (Et +0 ) (10) 
where 

o = tan! (A/B) (11) 
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Similarly, the fundamental component of the R.F. voltage may be 


written as 


where 


] 


a= atane 4 (C/D!) (13) 


Letting the electronic impedance es) be defined as 


Zo V(t) ¢/Ip(t) ¢ (14) 
then, upon using the relations given in Eqs. (10) to (13), the real 


and imaginary parts of Zz, are, 
Real (z,) = | V(t)¢/T,(t)F | cos($-c) (15) 


Imaginary ee) = | V(t)¢/T7(+)F | sin($-c) (16) 


In the computer program, Fourier analysis was carried out 
using over one hundred points for each period analyzed. As a check 
on the accuracy of the results, the waveform for the R.F. voltage 
was analyzed for the d.c. component and the harmonic components along 
with the fundamental. As the relation for V(t) was a cosine series, 
the analysis was considered accurate as long as the magnitudes of the 
d.c. and harmonic components were very smal] with respect to the 


magnitude of the fundamental component. 
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APPENDIX B 
DERIVATION OF REFLECTION COEFFICIENT FROM CARRIER REINSERTION 
The relationships between Voltage Standing Wave Ratio (VSWR), 
voltage reflection coefficient (IT), transmission line characteristic 
impedance hes and the load impedance (Z) are the following: 


ef 


SNR ee ee ep (1) 
(@) 


ear ee 
The magnitude of the voltage reflection coefficient is |r| where 


Iv. 


| 
LA 
and Vie V; are the reflected and incident voltage waves respectively. 


A further useful relationship is 


lv. + V| ke WOES 
elo IE eee es Sp (ie, (3) 
V5 Ny 


Defining the complex reflection coefficient as 


br fee ap |e) 
then 

Ir] =Ya2 +B? [rio = ae + Be (4) 
and 

+r | = (ita)? + Be 5 |ler[® = (14A)> + Be (5) 
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and 


Bs aitrtomsae (7) 


Once the values of +r | and ute are known, the values of A and B 

can be calculated. Both the magnitude and sign of A can be calculated 
however, only the magnitude of B can be calculated. The sign of B is 
arrived at by knowing the sign and magnitude of B at some time and then 
monitoring the value of fier]? For B to change sign the value of 


2 Sat ; 
|1+r|“ has to decrease to a minimum and then increase again. 


Both the relations jr4r|? and ir]? are a direct measure of 
the values of Iv.|4 and Iv If and therefore a direct measure of the 
incident and reflected powers. An R.F. detector operating in the 
square law region will give, as an output, a d.c. voltage that is 
directly proportional to the R.F. power incident upon it. Ideally then, 


an R.F. detector can be used to obtain the values of hier | and ir? 


as follows: 
MG ae 2 VitVnl§ 2 

| ~- |? - = |r+e|4 (8) 

Vila 

2 
Na 2 _ Vlg Ss 2 9 
Po i = TH = Irlg (9) 

IVilg 


where the subscript "d" refers to a quantity an R.F. detector can 
measure. Since the R.F. detector is not ideal, some proportionality 
constants will be required and the relations for the complex reflection 


coefficient will be 
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A=[ ujler[S - (148 [rio )1/2 (10) 
and 
B= as [ris - A (11) 


where U and S are arbitrary constants to be determined for each specific 


R.F. detector. 


In practice, the R.F. detectors used have characteristics 
that display appreciable variation when monitored over a wide range 
of incident R.F. signals. To account for a possible change in the 
values of the constants U and S, in Eqs. (10) and (11), the R.F. input/ 
d.c. output characteristics of the detector used, were measured over 
the range of R.F. signal levels of interest. These characteristics 
are shown in Fig. B.1 for various detector bias voltages. By making 
a piecewise linear approximation to these curves, values of U and S 


can be obtained for each individual linear approximation segment. 


The constants, U and S, are determined by measuring a load 
whose complex reflection coefficient is known. The general procedure 
is to take the d.c. output level of the R.F. detector in millivolts and 
find the equivalent R.F. power in milliwatts from Fig. B.1. This 
equivalent R.F. power is then the subscripted value ("d") in Eqs. (10) 
and (11). The values of U and S are those for the appropiate piecewise 
linear approximation that applies to that particular range of R.F. 
detected levels. This procedure is followed for the range of R.F. levels 
of interest. The complex reflection coefficient of an unknown load 
can be found by following a similar procedure of converting the 
detected d.c. levels to equivalent R.F. powers, from Fig. B.1, and 


substituting these values into Eqs. (10) and (11). 
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APPENDIX C 


R.F. DETECTION 
The R.F. detector unit used in the work reported here was 
a commercially available Hewlett-Packard Crystal Detector Model J424A. 
The waveguide mount has a bandwidth of 5.3 - 8.2 GHz, with the 
detector diode itself having a much broader frequency range. Schematic 
diagrams of the physical and electrical properties of this R.F. 


detector unit are shown in Fig. C.1. 


The crystal diode may be considered as a voltage generator 
with an internal resistance of 3K to 20K ohms shunted by approximately 
10 picofarads of capacitance. This gives the detector an R.C. time 
constant of from 30 to 200 nanoseconds. The equivalent internal 
resistance can be reduced to around 200 ohms by forward biasing the 
crystal diode; this places the value of the R C time constant in the 
range of 2 nanoseconds. Thus, the response time of the R.F. detector 
unit can be controlled by adjusting the forward bias of the crystal 


diode. 


This forward bias was applied in series with a 10K ohm 
resistor as shown in Fig. 6.5. As only a few hundred millivolts are 
required to forward bias the detector diode, a 10K ohm resistor was 
placed in series with the detector and the biasing d.c. power supply. 
This allowed the use of a commonly available power supply and also 
protected the delicate crystal diode from transients in the operation 
of the biasing power supply. The series 10K resistor was placed in 


the center conductor of a RG58 A/U coaxial cable and located inside 
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the BNC connector at the detector end of the cable. In this way 
parasitic effects were reduced while still providing a convenient 


method of applying the bias voltage to the detector unit. 


To ensure adequate response time in the detector system, 
the output of the detector was connected directly to the 50 ohm 
oscilloscope channel. This constituted a heavy loading of the R.F. 
detector unit which decreased the sensitivity of the unit and 
necessitated the use of the most sensitive range of the oscilloscope 
channel (10 mv/cm). By connecting the positive biasing voltage and 
the oscilloscope channel to the R.F. detector unit with a BNC tee, 
two advantages were realized. One advantage was that the actual for- 
ward bias across the detector could be monitored on the oscilloscope 
and also, this positive bias voltage offset the negative d.c. level 
that is the normal output of the detector unit, thus making it easier 


to use the most sensitive scale on the oscilloscope channel. 
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